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Partie A : Vibrations 

I - Vibrations libres des systèmes mécaniques à un degré de liberté 

1 - Système mécanique élémentaire 

2 - Amortissement visqueux 

3 - Méthode énergétique de Rayleigh 

4 - Notion de stabilité 

II - Réponse forcée des systèmes mécaniques à un degré de liberté 

1 - Réponse à une excitation harmonique 

2 - Différents types d'amortissement 

3 - Réponse à une excitation quelconque 

III - Systèmes à plusieurs degrés de liberté 

1 - Système à deux degrés de liberté 

2 - Généralisation aux systèmes à plusieurs degrés de liberté 

3 - Absorbeur dynamique 

4 - Introduction à l'analyse modale expérimentale 

Partie B : Acoustique 

IV - L'équation d'onde acoustique et ses solutions 

1 - Equation d'onde à une dimension 

2 - Equation d'onde en 3D 

3 - Ondes propagatives et ondes stationnaires 

4 - Intensité et puissance acoustique 

V - Notions d'acoustique physiologique et d'acoustique des salles 

1 - Introduction à la physioacoustique 

2 - Applications de l’acoustique statistique des salles 

VI - Mesures acoustiques 

1 - Chaine de mesure acoustique 

2 - Mesure des caractéristiques des matériaux 

3 - Mesure de la transparence acoustique 

4 - Mesure de la puissance acoustique 
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Documentation sur l'acoustique et les vibrations 



L'adresse des sites de l'Internet est indiquée, mais les changements d'hébergement sont 
encore fréquents. 

Publications d'organismes sur l'acoustique industrielle 

Brüel & Kjaer [http://www.bk.dk] 

Documentation sur les méthodes de mesure et d'analyse en acoustique et vibrations 
Acoustic noise measurements (5 th ed.) 

Noise control 

Frequency analysis (3 rd ed.) 

Brüel & Kjaer Technical Review (publiée régulièrement de 1949 à 1998 puis épisodiquement) 



Centre Technique des Industries Mécaniques (CETIM) [http://www.cetim.fr] 

Recueils de conférences et guides pratiques : 

Applications of Active Control to the Réduction of Noise and Vibrations (1996) 

Les nouveaux outils de la qualité acoustique. Comment façonner l'image sonore de vos 
produits (1996) 

Le bon usage du silencieux pour la réduction du bruit : machines, installations, véhicules 
(1997) 

Guide acoustique des installations de pompage (1997) 

Matériaux acoustiques pour l'industrie (2003) 

Réduction du bruit au voisinage des usines (2003) 

Centre d'information et de Documentation sur le Bruit (CIDB) [http://www.infobruit.org] 

° Annuaire Bruit des acteurs de l'environnement sonore édité chaque année (organismes 
publics, associations professionnelles, formations, laboratoires, organismes de contrôle, 
bureaux d'études, fabricants de matériaux acoustiques et d'instruments) 

° Echo-Bruit (magazine trimestriel de l'environnement sonore pour le grand public et les 
collectivités locales) 

° Acoustique & Techniques (revue trimestrielle destinée aux techniciens, acousticiens et 
tous professionnels concernés par les techniques de lutte contre le bruit) 



Centre Scientifique et Technique du Bâtiment (CSTB) [http://www.cstb.fr] 

° Cahiers du CSTB (revue technique dans le domaine du bâtiment avec des articles sur le 
bruit) 



Institut National de Recherche et de Sécurité (INRS) [http://www.inrs.fr] 

° Cahiers de notes documentaires - Hygiène et sécurité du travail (revue trimestrielle 
scientifique et technique avec des articles sur la prévention du bruit) 

Techniques de l'Ingénieur [http:/www.techniques-ingenieur.fr], voir en particulier la rubrique 
"mesures et contrôle" 



Revues internationales 

Les sites web permettent la plupart du temps d'avoir accès aux résumés des articles des 
dernières parutions. 

Applied Acoustics [http://www.elsevier.eom/inca/publications/store/4/0/5/8/9/0/] 

Acustica United with acta acustica [http://www.icp.inpg.fr/ACTA/] 

International Journal of Acoustics and Vibration [http://www.rcom.ru/IJAV/] 

Journal of Acoustical Society of America [http://ojps.aip.org/jasa/] 
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Journal of Sound and Vibration [http://www.academicpress.com/jsv/] 

Journal of Vibration and Acoustics (Transactions of ASME) 
[http://www.asme.org/pubs/journals/vibrate/vibrate.html] 

Noise Control Engineering Journal [http://users.aol.com/inceusa/indxncej.html] 

Principales sociétés d'acoustique 

Société Française d'Acoustique (SFA) [http://www.loa.espci.fr/sfa/] 

Association Française de Mécanique (SFM) [http://www.afm.asso.fr/] 

European Acoustics Association (EAA) [http://eaa.essex.ac.uk/eaa/] 

International Institute of Noise Control Engineering (IINCE) [http://users.aol.com/iince1/] (et INCE 
USA [http://users.aol.com/inceusa/]) 

International Institute of Acoustics and Vibration (IIAV) [http://www.iiav.org] 

Acoustical Society of America (ASA) [http://asa.aip.org] 



Diverses sources de documentation et d'information 

Association Française de Normalisation (AFNOR) : normalisation, certification, catalogues de 
normes, directives, formations, forums [http://www.afnor.fr] 

Centre d'information et de Documentation sur le Bruit (CIDB) : publications, bibliothèque 
documentaire, stages, colloques. Edite chaque année l'Annuaire des acteurs de l'environnement 
sonore (organismes publics, associations professionnelles, formations, laboratoires, organismes de 
contrôle, bureaux d'études, fabricants) [http://www.cidb.org] 

Laboratoire National d'Essais (LNE) métrologie acoustique, étalonnages, certification 
[http://www.lne.fr] 

Institut de l'Information Scientifique et Technique (INIST-CNFtS) : fourniture de documents 
scientifiques et techniques [http://www.inist.fr] 



Réglementation 

Les Journaux Officiels [http://djo.journal-officiel.gouv.fr/] chercher avec le mot clé "Bruit". En 
particulier : 



° Législation communautaire en matière d'environnement /Volume 5 - Bruit 
Situation au 01.09.91 (Edition de 1993 - réf. : 3699100541) 

Situation au 30.06.94 (Edition de 1996, complément à I' Edition de 1993 - réf. : 
3699600541) 



° Bruit - Prévention, maîtrise et contrôle des nuisances sonores (réf. : 31 3830000, édition du 
2 juin 1995) 

Cet ouvrage regroupe l'ensemble des textes législatifs et réglementaires 
d'origine nationale et communautaire qui régissent la prévention, la maîtrise et le contrôle 
des nuisances sonores. 
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I - VIBRATIONS LIBRES DES SYSTEMES MECANIQUES 
A UN DEGRE DE LIBERTE 



1. SYSTEME MECANIQUE ELEMENTAIRE 

1.1 - Equation du mouvement 

1.2 - Raideurs équivalentes 

1.3 - Masses équivalentes 

1.4 - Linéarisation 

Application 1 : Détermination de la fréquence propre 

2. AMORTISSEMENT VISQUEUX 

2.1 - Mouvement sous-amorti 

2.2 - Mouvement sur-amorti 

2.3 - Mouvement avec amortissement critique 

2.4 - Equation normalisée d’un système à 1 degré de liberté 
Application 2 : Mesure du taux d'amortissement 1 



3. METHODE ENERGETIQUE DE RAYLEIGH 

3.1 - Equation du mouvement 

3.2 - Pulsations propres 

3.3 - Exemples 



4. NOTION DE STABILITE 



ANNEXE : Equations différentielles et mouvement harmonique 
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I - Vibrations libres à 1 ddl 



1 - SYSTEME MECANIQUE ELEMENTAIRE 

Le modèle du système mécanique élémentaire considère le mouvement d'une masse (corps 
rigide) par rapport à une partie fixe. On parle de système à 1 degré de liberté (DDL) 
quand la masse unique a un mouvement dans une seule direction (translation ou rotation 
autour d'un axe). Si une ou plusieurs masses ont des mouvements de translation et de 
rotation dans plusieurs directions, il s'agit d'un système à plusieurs degrés de liberté. 
Malgré sa simplicité le système à 1 DDL peut représenter le comportement dynamique 
de systèmes très variés dans le domaine des basses fréquences. La modélisation considère 
une masse équivalente en mouvement qui possède des liaisons avec les parties fixes 
caractérisées par une raideur équivalente (souvent schématisé par un ressort). 



1.1 - Equation du mouvement 

Il existe principalement deux types de systèmes : les systèmes de translation et les 
systèmes de torsion. 



1.1.1 - Systèmes avec mouvement de translation 

Ils sont schématisés par le système masse - ressort : la masse m [en kg] est animée d'un 
mouvement de translation dans la direction x auquel s'oppose la force due à la raideur du 
ressort. Dans le domaine linéaire du ressort, le coefficient de raideur k [en N/m] est une 
constante et la force de réaction F k = —kx . 



Le principe d'Alembert permet d'écrire l'équilibre dynamique du système masse-ressort 
(. Figure 1.1) 

F k (t) = F i (t) 



entre 

et 



la force dynamique F i (t) = m — : 

dr 

la force de réaction exercée par le ressort F k (/) = —kx 
(le signe - est du à la force qui s'oppose au déplacement) 



= mx 



-> X 



WW m 





F = mx 



Figure 1.1 - Equilibre du système masse-ressort 



L'équation différentielle du mouvement (homogène du second ordre) se déduit de 
l'équation d'équilibre précédente 



mx + kx = 0 



( 1 ) 
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Sa solution (voir Annexe) est une fonction périodique pour le déplacement 

x(t) = A sin(&V + </>) (2a) 

A et (f) sont l'amplitude et la phase qui dépendent des conditions initiales et co est la 
pulsation (ou fréquence angulaire) 



vitesse : 



accélération : 



dx(t) 

dt 

d 2 x(t) 

dt 2 



= x = coA cos (cot + (/)) 



= x = -coAsm(cot + (j)) 



(2b) 

(2c) 



L'équation (1) s'écrit alors 



-mû) 2 A sin(&V + </>) = —kA sin(&V + (fi) 



La pulsation naturelle est la valeur de a qui satisfait la relation 




[en rad/s] 



( 3 ) 



û) 0 ne dépend que des constantes mécaniques du système. Les constantes A et (j) 
dépendent des conditions initiales (x est intégré deux fois pour obtenir x, d'où deux 
constantes d'intégration). Les conditions initiales sont le déplacement x 0 et la vitesse v 0 à 
l'instant t = 0 : 



x 0 = x(o) = Asin(ft> 0 0 + (/)) = Asin^ [m] 

v 0 = x(o) = ty 0 Acos(&> 0 0 + ^) = <y 0 Acos^ [m/s] 



donc û)qxI + vf = a>l A 2 (sin 2 (f) + cos 2 (f>) et S ' n ^ 

cos^ 



Û) 0 X 0 



, soit 



A = 






,.2 2 . 2 
Û) Q x 0 + v 0 

û) n 



et </> = arctan 



û) Q x 0 



( 4 ) 



Le déplacement 




( 5 ) 



est la réponse libre du système à 1 degré de liberté non-amorti. 

C'est une fonction harmonique à la pulsation naturelle û) 0 dont l’amplitude est imposée par 
les conditions initiales. Cette amplitude reste constante car la modélisation n'a pas pris en 
compte le phénomène de dissipation d'énergie présent dans tout système mécanique. Dans 
la réalité, une décroissance de l'amplitude avec le temps sera observée. 
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I - Vibrations libres à 1 ddl 



T = Itt/cOq 





Figure 1.2 - Réponse libre du système à un degré de liberté ( déplacement et vitesse) 
de pulsation propre (O, = 2jtf {) . 



Remarque 1 

La solution de l'équation différentielle peut aussi s'exprimer par la somme d'une 
fonction sinus et d'une fonction cosinus : 

x(t) = A sin(ft> 0 t + (/)) = A sin (/) cos co Q t + A cos (j) sin co,t = B l cos û) 0 t + B 2 sin û) 0 t 

En utilisant cette forme de solution, les conditions initiales conduisent à 

v n 

x{t) = x, cos O) 0 t H sin co Q t . 

co. 

Il est aussi possible d'exprimer la solution à l'aide d'une seule fonction cosinus 

x{t) = Acos(co 0 t + 0) = Acos 0 cos co 0 t + Asin#sinty 0 t 

A partir de l'expression précédente 

Acosé? = x 0 

A siné? = - v 0 /ft» 0 

ce qui permet de retrouver la même amplitude que précédemment 
A = sJ A 2 cos 2 6 + A 2 sin 2 6 = -^x\ +v\f(ù\ 

et conduit à la phase 

„ - v, 

0 = arctan 

co 0 x o 

On note qu'il existe une relation entre 0 et <j> l'angle de phase associée à la solution 
sinus car 

tan 0 = — = - cot (j) d'où 0 = (j) - n . 

co 0 x 0 
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Remarque 2 




t F k=~ k ( x o+*) 



T F t = mx + mg 



Figure 1.3 - Système masse-ressort suspendu 



Système suspendu : au repos, la masse exerce une force de traction mg sur le ressort qui se 
traduit par un étirement initial x 0 . L'équation du mouvement s'obtient à partir de l'équilibre 
dynamique des forces appliquées à la masse : F k ( t ) - F j ( t ) = 0 

mx + mg + kx + kx 0 = 0 

Avec mg = -kx 0 , on obtient la même équation du mouvement x(t) que précédemment 
(centrée sur la position d'équilibre x 0 ). 



1. 1.2 - Systèmes de torsion 

Un corps rigide oscille autour d'un axe (vibration de torsion). 




Figure 1.4 — Equilibre d'un système de torsion. 



L'équilibre 
de torsion 



dynamique entre le moment dynamique M t ( t ) = 7, 



d 0 

dt 2 



= Jn0 



M t (t) = —k t 6 permet d'écrire l'équation du déplacement angulaire 



et le moment 

d{t) 



J q0 + k t 0 — 0 



J 0 moment d'inertie de la masse et k t raideur de torsion [en Nm/rad]. 
La pulsation naturelle est 



co 



0 




[rad/s] 



( 6 ) 



(7) 
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I - Vibrations libres à 1 ddl 



et la réponse libre du système non-amorti 



9{î)-. 



) 2 

7 0 



sin 



û) n 



O) 0 t + arctan 



Exemple : pour un disque de diamètre D et 
d’épaisseur h , à l’extrémité d'un 
arbre de longueur 1 et de diamètre d 



Jn 



mD 1 
~8~ 



et k, 



G7ld- 



32/ 



avec la masse m = ph 



kD 



et le 



module de cisaillement G. 



&o j 




(B) 



A 

V 



1.2 - Raideurs équivalentes 

Le modèle le plus simple pour représenter la raideur correspond au ressort suspendant une 
masse. En négligeant la masse du ressort, les forces agissant sur la masse sont 

■ la force de gravité F = mg 

■ la force de réaction du ressort F k = -kx 

Le travail effectué en déformant le ressort est stocké dans le système sous la forme d'une 
énergie potentielle (travail de la force: F - kx) 



U = jka da 

o 




Le coefficient de raideur k est une constante qui montre que la force de réaction est 
proportionnelle à la déformation (qui augmente avec la masse). Ce modèle est valide 
jusqu'à un certain point au-delà duquel k n’est plus constant. 



k / repos 

FWW 



x 

0H>: 



IhA/VV 



-> F 




Figure 1.5 - Zone de comportement linéaire d'un ressort de constante k = 750 N/m. 
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Figure 1.6 - Trois configurations pour un système masse ressort 



1.2. 1 - Raideurs en parallèle 




F F F 



Figure 1.7 -Raideurs en parallèle 
F = k\X + k 2 x = k eq x 

La raideur équivalente est la somme des raideurs k 1 et k 2 

k eq =^i+^2 

Ce résultat se généralise pour n raideurs en parallèle 

k eq = k l +k 2 +... + k n (10) 
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I - Vibrations libres à 1 ddl 



1.2.2 - Raideurs en série 




F F 



Figure 1.8 - Raideurs en série 



L'équilibre des forces conduit aux relations suivantes : 

7 7 k eq 

k.x. = kx => X. = X 

l I et/ l 

7 7 

k 2 X 2= k ec X => X 2=~j— X 

k 2 



L’élongation totale s’écrit alors 
d'où 




Ce résultat se généralise pour n raideurs en série 



111 1 
— = — + — + ... + — 



k l k 2 



1.2.3 - Raideurs de flexion 

Dans l’ouvrage de Rao (voir Références), des raideurs de flexion (translation) et de torsion 
sont données pour quelques systèmes mécaniques simples. Elles peuvent être obtenues en 
utilisant les résultats du cours de RDM. Par exemple, la déformation statique d'une poutre 
sur appuis simples chargée par une masse M se calcule par 




Figure 1.9 - Raideur équivalente d'un système masse-poutre 
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La raideur équivalente de la poutre s'obtient par ( k eq S = Mg = P ) 



où 8 = y(o) 
d'où 



Pa 2 b 2 
3 EIL 



k 



eq 



P _ Mg 

~S~~S~ 



est la déflexion de la poutre à l'emplacement de la masse ( x = a ), 



k 



eq 



P 

~s 



3E1L 

a 2 b 2 



1.3 - Masses équivalentes 



Principe : Si T correspond à l'énergie cinétique totale de toutes les masses en mouvement 
(en translation et en rotation), 

T = E;Ui, 2 +XL ,è) 

i Z j Z 

la masse équivalente est définie par 




x est la vitesse au point où on veut calculer la masse équivalente (généralement, où se 
trouve définit la raideur équivalente). 

Note : 

a) les masses en mouvement doivent être connectées rigidement entre elles sinon il 
s'agit d'un système à plusieurs degrés de liberté. 

b) pour les systèmes en mouvements de rotation, on considère le moment d'inertie 
équivalent 

1 

t = -j e 2 
2 ^ 



1.4 - Linéarisation 
1.4. 1 - Pendule simple 



j 0 ë=M t 

M t = F t l = -mg l sin 0 



Figure 1.10 - Pendule simple 
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I - Vibrations libres à 1 ddl 



Equation différentielle : J o 0 + mgl sin 0 = 0 

Linéarisation pour les petits mouvements: sin 0-0 => J {j 0 + mgl 0 = 0 

avec J 0 =ml 2 

ml 2 0 + mgl 6 = 0 

La pulsation propre du système est 





G : centre de gravité 

J Jd + Mgcl sin# = 0 

J Jd + Mgd 0 = 0 (Eq. linéarisée) 




Figure 1.11 — Pendule composé 



j n 

Longueur équivalente du pendule simple : / = 

Md 



[APPLICATION 1] Détermination de la fréquence propre 

Détermination de la fréquence propre à partir de l'écrasement x 0 du système de suspension 




Figure 1.12 - Fréquence propre des plots anti-vibratiles 
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Vibrations - Acoustique 1 



Ce principe s'applique également aux systèmes continus : ici une masse fixée sur une 
poutre. Si la masse de la poutre est négligeable devant la masse de l'ensemble, la poutre 
n'apportera que de la raideur au système (mouvement de flexion). La déflexion statique S 
observée à l'emplacement de la masse permet d'obtenir la première fréquence propre du 
système 




Par exemple, pour le système de la Figure 1.9 : C0 Q 




2 - AMORTISSEMENT VISQUEUX 




Figure 1.13 — Exemple d'amortisseur 

L'amortisseur dissipe l'énergie. L'écoulement laminaire de l'huile crée une force de réaction 
proportionnelle à la vitesse du piston 

F _ Mt) 

c dt 

Le coefficient d'amortissement c dépend de la viscosité [en Ns/m] 

Autre type d'amortisseur sensiblement proportionnel à la vitesse : un bloc de caoutchouc. 




Equilibre des forces : 

F k + F c - F, = 0 

d'où 
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I - Vibrations libres à 1 ddl 



l'équation décrivant le mouvement d'un système amorti à 1 degré de liberté. En posant 
x(t) = a e" , cette équation devient (mr 2 + cr + k)a e rt =0. Puisque a e" ne peut pas être 
nul quel que soit t , c'est l'équation caractéristique suivante qui doit être vérifiée 

mr 2 + cr + k = 0 . 

Ses solutions sont 

r\ j _ c a/c 2 - 4 km 

r 2 J 2m 2m 

Plusieurs types de solutions sont envisageables en fonction de la valeur du discriminant 
c 2 - 4 km : deux racines réelles, une racine double ou deux racines complexes. 

Pour faciliter l'analyse, on définit le coefficient d'amortissement critique c cr 

c 2 - 4 km = 0 => c cr = 2yfkm = 2 mû) 0 

co 0 pulsation naturelle non-amortie (fréquence ou pulsation propre) 

Définition du taux d'amortissement 

Ç _ c _ c 
c cr 2 mû) 0 

Les racines de l'équation caractéristiques peuvent donc s'écrire en fonction du taux 
d'amortissement 



1 1 = -Ç(0 Q ± co 0 Jç 2 -1 x(t) = a x e' 1 ’ + a, e r - 

r 2 J 

En fonction de la valeur de Ç , trois types de mouvement peuvent être observés : 
le mouvement sous-amorti (oscillations vibratoires amorties), 
le mouvement sur-amorti (retour à la position d'équilibre sans oscillations), 
le mouvement avec amortissement critique, qui correspond à la limite entre les 
deux cas précédents. 

2.1 - Mouvement sous-amorti ( 0 < Ç <1 ) 



n =-£»<> -Ma/ï-T 



= -Cû ) 0 + jcoJl-Ç 2 



avec j 



r = -v^r et vr-i=V(-i)(i-<r 2 )=-À/i-r . 



La solution de l'équation différentielle est de la forme 



x 



(t) = 



+a 2 e~ j ^ 

v 






a x , a 2 sont des constantes complexes arbitraires 



x(t) = A e - û>0 ' sin(ty/ + (j>) 
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où A et (j) sont les constantes d'intégration et cû d est la pseudo-pulsation (fréquence 
naturelle amortie) 




x 0 O) d 

</> = arctan — . 

v 0 + Ça> 0 x 0 




temps [s] 

Figure 1.14 - Réponse libre d’un système sous-amorti ( Ç < 1 ) 

2.2 - Mouvement sur-amorti ( Ç >1 ) 

h = ~Çg>o + «o yjC 2 -1 F = -Ç® o - JC 2 -1 

La solution pour le déplacement est 




a j et a 2 sont réels. 

Les conditions initiales 



x 0 = x(û) = a 1 + a 2 
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Un système sur-amorti n’est pas un oscillateur 



I - Vibrations libres à 1 ddl 






temps [s] 



Figure 1.15 - Réponse libre d’un système sur-amorti ( Ç > 1 ) 
a) x 0 = -1 et v 0 = 0 , b) x 0 = -0.5 et v 0 = 3 , c) x 0 = 0 et v 0 = 5 . 



2.3 - Mouvement avec amortissement critique (^=1 ) 

C’est la valeur de Ç qui sépare le mouvement oscillant d’un mouvement non-oscillant. La 
racine double fournit une solution particulière 

'i = h = = ~a> 0 ■ 

La solution générale de toute équation différentielle du second ordre est donnée par la 
combinaison linéaire de deux solutions particulières indépendantes. La première solution 
est 

-(Ont 

x l =e 0 

L’intégration complète de l’équation différentielle se fait en définissant une deuxième 
solution de la forme x 2 = z(t) e~ aof =(a+ fit) e~ w ° l , ce qui conduit à la solution générale 

x(t) = x 1 + x 2 = (a l + a 2 t) e~ a ° t 
Les conditions initiales permettent d’identifier les constantes 
x 0 = a 1 

v 0 = (- û) 0 a l + a 2 e~ aot - û) Q a 2 t e~ c ° ot ) = ~( 0 Q a l + a 2 
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est d’obtenir la solution cherchée 




temps [s] 



Figure 1.16 - Réponse libre d’un système ayant un amortissement critique ( Ç = 1 ) 
a) x 0 = 0.5 et v 0 = -4 , b) x 0 = 0.5 et v 0 = 4 . 



2.4 - Equation normalisée d’un système à 1 degré de liberté 

L’équation différentielle du mouvement d’un système à 1 ddl est souvent représentée avec 
des termes normalisés par rapport à la masse m : 

c . k 

x H x H x = 0 

m m 

]ç Q Q 

En considérant — = &> ( , et puisque Ç = , on obtient la relation — = 2û) 0 Ç qui 

m 2mû) 0 m 

conduit à la forme normalisée de l’équation différentielle 

x + 2Ç û) 0 x + (Oq x = 0 



[APPLICATION 2] Mesure du taux d’amortissement 

Le coefficient d’amortissement ou le taux d’amortissement sont les plus délicats à 
déterminer. Si la masse et la raideur peuvent être déterminées par des tests statistiques, 
l’amortissement nécessite une mesure dynamique. 

Une approche consiste à mesurer la décroissance de l’enveloppe pour un système sous- 
amorti : les points de mesure x(t { ) et x(t 2 ) correspondent à A <? _Wl et à A . 
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Cette approche conduit au concept de décrément logarithmique : 



8 = ln 



x{t) 

x(t + T) 



^ 2.71 

avec T - — 
CO , 



la pseudo-période d’oscillation. 




Figure 1.17 — Réponse libre amortie 



avec x(t) = A e w sin (co d t + <f>) . Puisque co d T = ln , 

x(t + T)= A e~^ a>o{ ' ,+T ' 1 sin (co d t + û) d T + </>)= A e ~^ 0>0 ' e ~^ T sin(ty^t + (/)) 

on obtient 

8 = ln e C ° ),,r = Çco 0 T . 



2 n 

En posant T = = 



co. 



, le décrément logarithmique s’écrit S = et 






CO, 






permet d’exprimer le taux d’amortissement 




Le décrément est obtenu par des mesures en t ] et t 2 par 8 = ln — — - 

x[t 2 ) 

Ç il est aussi possible de calculer le coefficient d’amortissement 



. Connaissant m , k et 



c = 
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3 - METHODE ENERGETIQUE DE RAYLEIGH 

x(t) = A sin(&J 0 f + </>) 
x(t) = co Q Acos(co 0 t + </>) 

Figure 1.18 - Déplacement et vitesse d'un système à 1 ddl en vibration libre 

En faisant l'hypothèse qu'il n’y a pas d’énergie dissipée : 

Energie potentielle en fonction du temps 

X J J 

U{t)= \ka da =—k x 2 = —k A 2 sin 2 (co 0 t + (/)) 
o 2 2 

Energie cinétique en fonction du temps 

T(t) = — mx 2 = — m cüq A 2 cos 2 (<» 0 t + ^) 

1 . 2 

En considérant k = mco ( j , 

U{t) + T{t ) = ^mû)l A 2 [sin 2 (û) 0 t + (/)) + cos 2 ( a> Q t + ç. i)] 

1 ,.2 a 2 1 , ,2 

= —m(OnA =—kA 

2 2 

L’énergie totale instantanée du système conservatif en mouvement est proportionnelle à 
A” et indépendante du temps t : 

U{t) + T{t) = Cte, Vr 

Deux conséquences importantes : 

1) La dérivée de l’énergie totale (potentielle + cinétique) est nulle, d'où la relation RI 




2) A deux instants t l et t 2 quelconques 

U{t l )+T{t l ) = U{t 2 )+T{t 2 ) 
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1 - Vibrations libres à 1 ddl 


En choisissant t ] et t 2 


, tels que 






si 


U{h) = ^max 


alors 


T{t i) = 0 


si 


U{t 2 ) = 0 


alors 


) = ^max 



on obtient la relation R2 




valable pour les systèmes conservatifs soumis à des mouvements harmoniques. Ce 
principe de conservation est utilisé : 

• pour obtenir des équations du mouvement, 

• pour calculer directement la pulsation naturelle des systèmes. 

3.1 - Equation du mouvement 




Figure 1.19 - Masse suspendue en vibration libre 



L’énergie potentielle totale du système est la somme de l’énergie potentielle élastique du 
ressort et de l’énergie potentielle gravitationnelle. Toutes les deux doivent être exprimées 
comme un changement à partir d’une position arbitraire. On choisit ici (voir Figure 1.19 ) 
la position d’équilibre x 0 qui correspond à une énergie élastique stockée dans le ressort 

U t =lkx o 2 . 



L’énergie potentielle U 2 à une distance x de la position d’équilibre correspond à une 
augmentation de l’énergie élastique qui devient -^-k(x 0 +x) 2 et une diminution de 
l’énergie potentielle mgx due au changement de position 



U 2 = ~^k(x 0 + xf - mgx . 



La différence d’énergie potentielle est 



1 2 1 1 

U = U 2 - U l = —I<.Xq + —kx~ + k x 0 x - mg x - — k x 



2 

0 
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et en considérant que kx 0 = mg , il reste 



U = 




L’énergie totale (cinétique + potentielle) est 



1 2 1 ^ 

T + U = — m x H — k x = Constante 

2 2 



et sa dérivée par rapport au temps (RI) 

d(T + U) 
dî 



m xx + k xx = 0 



Puisque la vitesse x n’est pas identiquement nulle (quel que soit t) dans l’expression 
précédente, on peut en déduire l’équation différentielle du mouvement 



m x + k x = 0 



3.2 - Pulsation propre 

La méthode énergétique peut aussi être utilisée pour obtenir la fréquence propre (naturelle) 
des systèmes oscillants conservatifs. Il suffit d’appliquer la relation : 



U =T 

max max 



1 



— kxZ 



— Tî IV 

2 max 2 * v ■' v max 



Pour un mouvement périodique d’amplitude A et de pulsation propre co 0 , x max = A et 
v max =« 0 A,d’où 



ce qui conduit à 



-mû ) q A 2 = —kA 2 
2 0 2 



û) 0 




3.3 - Exemples 
3.3. 1 - Pendule simple 

La masse est supposée quasiment ponctuelle par rapport à la longueur / du bras et la 
masse du bras est considérée comme négligeable. Dans ces conditions, le moment d’inertie 
autour de l’axe est 

J 0 = ml 2 



20 




Vibrations - Acoustique 1 



I - Vibrations libres à 1 ddl 



8 

Figure 1.20 - Pendule simple 



Le déplacement angulaire est mesuré à partir de la position d’équilibre. L’énergie cinétique 
du système est 

T = -J o 0 2 = —ml 2 9 2 , 

2 0 2 

et l’énergie potentielle 

U = m g (l - h) = m g l{ 1 - cos 6 ) . 

La dérivée de l’énergie totale s’écrit alors 

d{T + U) d 1 .i ■-> . 

= — —ml 0 + m g l (1- cosÛ) =0 

dt dt\_ 2 

ml 2 66 + m g l (sin 0)6 = 0 
6 (l 6 + g sin 0 ) = 0 . 

La vitesse angulaire d(t ) ne peut être nulle à tout instant donc 0 + y sin 0 = 0 et après 

linéarisation pour les petits mouvements (sin 0 ~ 0) 

6+^6=0 

l 

Cette équation permet de calculer le déplacement angulaire et la pulsation naturelle 
û ) o = sj g/l du pendule. 

Note : La méthode énergétique peut aussi s’appliquer directement sur 

0 ( 0 \^ 0 2 

l’approximation des petits déplacements : 1 - cos 9 = 2 sin 2 — = 2 — = — . 

Pour obtenir directement la pulsation propre par la méthode énergétique, il faut écrire 

2 

=lml 2 (w o e„) 2 U m „ =mglY 

tj —t m 2 ~ — 

u max — ' max ~ ^ 

3.3.2 - Ressort pesant 

Détermination de la pulsation naturelle d’un système masse-ressort dont la masse du 
ressort n’est pas négligeable devant celle de la masse. 
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L’énergie potentielle ne 
change pas , mais la masse 
du ressort m r doit entrer 
dans le calcul de l’énergie 
cinétique. 




4 - NOTION DE STABILITE DES SYSTEMES 

Dans ce qui précède, les coefficients m , c et k sont considérés comme positifs. Les 
solutions du mouvement peuvent se classer en 4 catégories : libre ( c = 0 ), sous amorti, sur- 
amorti et avec amortissement critique. L’amplitude des mouvements est finie. Le système 

est stable. 
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L’Annexe montre que dans le cas où m > 0 , c = 0 mais k < 0 , le déplacement solution du 
système est de la forme 

x(t) = A sinh û)t + B cosh cot . 

L’amplitude du mouvement s’accroît avec t . De telles solutions sont appelées divergentes 
ou instables. 




■ 



vwwiw/wf 



■ * 



V 



2 



/ sin 6 

l< > 



J 

<r 



2 ky 

< 



x 




Figure 1.22 - Exemple de système instable : le pendule inversé 



kl 

On a 2ky = 2 — sin 6 
2 



Energie cinétique : 



T = -jJ 2 =-ml 2 è 2 
2 ° 2 



Energie potentielle du ressort 



U r = — (2k)y 2 avec y = — sin# 
r 2 2 



U =- 
2 



f U 2 • 2 r\ 

sin 0 

2 



Energie potentielle de gravitation: 



avec x = l-l cos 0 



U g = —m g x 

= mgl (cos 0- l) 



L’énergie totale est 



1 , 1 

T + U = —ml 2 0 2 +- 
2 2 



(k{1 ' 

sin 0 

2 



+ m g l (cos 6 — l) = Constante 



et sa dérivée 



d 



2 

— {T + U) = ml 2 6Ô + — sin#cosé?#-mg/sin## = 0 
dt 2 



Ce qui conduit à 



t •• kl 

ml 6 + — — sin 6 cos 6 — m g l sin# = 0 
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et avec l’hypothèse des petits mouvements 6 «1 => sin# = # et cos#~l, à 
l’équation linéarisée du mouvement 



ml 0 + 




0 = 0 



critère de stabilité : — - mg > 0 . Dans le cas contraire, le système est instable et le 

mouvement diverge 
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ANNEXE 

Equations différentielles et mouvement harmonique 



La variable y est une fonction de z et solution de l’équation différentielle homogène du 
second ordre 



cl 2 y 

dz 1 



+ ay = 0 



En supposant que la solution est de la forme y = ae rz : 



r7 d y 
— = are , 



dy 

dz 



dz 



ar 2 e rz et ( r 2 +a)ae rz =0 



La solution non-triviale (y ^ 0) est donnée par les solutions de l’équation caractéristique 

r 2 + a = 0 

qui possède deux racines réelles ou imaginaires selon le signe de a 
1) a négatif : a = -co 2 
Les racines sont 
et la solution 



y" - co 2 y = 0 et r 2 - CO 2 = 0 



q = CO et r 2 = -CO 



y = Ae +Be ou y = C cosh co z + D sinh co z 



(car e az = coshttfz + sinhtüz et e coz = cosh coz - sinh coz ) 



'y 

2) a positif : a = (0~ 
Les racines sont 
et la solution 



y" + û) 2 y = 0 et r“ + &r = 0 
r { = jCO et r 2 = —j(û 



y = Ae jaz +Be~ jû,z . 



A et B sont des constantes complexes. Dans le cas où une solution réelle est recherchée, 
la partie imaginaire doit être nulle 

y = A (cos (O z + j sin co z) + B (cos co z - j sin co z) 

= (A S + jAj ){coscoz + jsincoz)+ {B r + jB r )(coscoz - jsincoz) 

= (A r +B r )cos coz~{A i - B j )sin co z + j [(A 7 + B 1 )cos co z + (A r - B R )sin co z] 



Les relations suivantes entre A et B annulent la partie imaginaire 

Aj = -B j et A r = B r , soit B = A* 

et permettent d’écrire la solution réelle sous les formes 

y = Ae J0>z + A* e~ ,az = 2 R e{Ae JÜ>z }= 2 A R cos coz — 2 A 1 sinty^ 
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.V 



= C cos a z + D sin co z ou y = E sin(ft) s + i 



car E sin(&> z + </>) = E sin ^ cos û)z + E cos <f> sin co z , donc 

C = E sin </> et D = E cos </) 
et 

E = Vc 2 + D 2 et à- arctan — . 

C 



Résumé 

Equation harmonique : ^ + co 1 u(t ) = 0 

dr 

Différentes formes équivalentes de la solution : 

u{t)=Ae iû> ' +Be- j0) ' avec B = A* et A = C/2-jD/2 
u(t) = C cosat + Dsinœt avec C = Esin<f> et D = Ecos(/> 

u(t) = Esin(cot + (f)) avec E = v/c 2 + D 2 et (f) = arctan — 
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Il - REPONSE FORCEE DES SYSTEMES MECANIQUES 
A UN DEGRE DE LIBERTE 



1 . REPONSE A UNE EXCITATION HARMONIQUE 

1.1 - Excitation d’un système non-amorti 

1.2 - Notation complexe 

1.3 - Excitation harmonique des systèmes amortis 

1 .4 - Excitation par la base 

1.5 - Excitation par déséquilibre dynamique en rotation 
Application 1 : Accéléromètre 

Application 2 : Isolation vibratoire - techniques de base 



2. DIFFERENTES FORMES D’AMORTISSEMENT 

2.1 - Amortissement de Coulomb 

2.2 - Energie dissipée par cycle et amortissement équivalent 

2.3 - Autres formes d'amortissement 



3. REPONSE A UNE EXCITATION QUELCONQUE 

3.1 - Réponse impulsionnelle 

3.2 - Réponse à une excitation non-sinusoïdale 

3.3 - Réponse à une excitation périodique quelconque 

3.4 - Utilisation de la transformation de Laplace 

3.5 - Réponses à des excitations aléatoires 
Application 3 : Analyseur de spectre 

Application 4 : Mesure des fonctions de transfert et de la 
réponse impulsionnelle 
Application 5 : Mesure de l'amortissement 2 

Annexe - Informations pratiques sur les mesures vibratoires 



1107 




I I - Vibrations forcées à 1 ddl 



Vibrations - Acoustique 1 



28 




Vibrations - Acoustique 1 



I I - Vibrations forcées à 1 ddl 



1 - REPONSE A UNE EXCITATION HARMONIQUE 

On étudie la réponse d’un système amorti à 1 ddl à une excitation harmonique sinusoïdale 
produite par une force extérieure au système. Ce type d’excitation se rencontre 
fréquemment dans l’industrie (machines tournantes, ventilateurs, moteurs, pompes ...). Les 
résultats obtenus avec une excitation sinusoïdale pourront s’étendre à des excitations 
harmoniques plus complexes représentées par des séries de Fourier, en appliquant le 
principe de superposition. On suppose donc qu’il s’agit de systèmes linéaires. 



1.1 - Excitation d’un système non-amorti 




c 



F(f) 



Equation du mouvement en x 



Dans un premier temps, l’amortissement 
est considéré comme négligeable (c = 0). 
La force d’excitation harmonique est 

F(t)= F cos cot 

F amplitude crête, Cû pulsation de 
l’excitation ou pulsation forcée. 



m x(r) + k x(t ) = F cos co î 
normalisée par rapport à la masse 

P 

x(t) + co: x(t) = — cos cot = f cos cot 
m 



Cette équation est une équation différentielle linéaire non-homogène et sa solution est la 
somme de la solution de l’équation homogène ( / = 0) et d’une solution particulière. La 
solution particulière peut souvent être obtenue en supposant qu’elle a même force que la 
fonction d’excitation 

x p = A cos cot . 

En l’introduisant dans l’équation du mouvement on obtient 

- co 2 A cos cot + col A cos cot = f cos cot 

d’où on tire 

f 

A = — — tant que coïû) 0 . 

col ~ (0 ~ 

En ajoutant cette solution particulière de l’équation non-homogène à la solution générale 
de l’équation homogène, la solution générale du système est 

f 

x(t) = A l sin (O 0 1 + A 7 cos co 0 t H — — cos cot . 

col ~ ûr 
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Les coefficients A l et A 0 sont déterminés à partir des conditions initiales 

/ 



x(0) = x 0 = A 2 + 

a>- - co 

x(0) = v 0 = û) 0 A! . 



La réponse forcée est donc 



c(?) = -^-sin<y 0 t + 



0)a 



f 



* 0 ~ 2 2 

V (°o ~C° J 



COS CO n t + 



/ 



o l 7 

O) o - O) 



cos cot . 



Le mouvement est très dépendant des conditions initiales, comme le montre la figure ci- 
dessous. 




temps [s] 

(a) 




temps [s] 

(b) 



Figure 2.1 — Réponse du système non-amorti à une excitation harmonique dont la 
pulsation est le double de sa pulsation naturelle : (a) conditions initiales x () et v 0 

non nulles, (b) condition initiale v 0 nulle. 



Quand la pulsation d'excitation co tend vers co 0 un phénomène de battement apparaît (à 
partir de l'expression précédente du déplacement, il est possible de montrer qu'il peut 
s'écrire comme le produit de deux fonctions harmoniques des pulsations (O + (O 0 et 

co - co Q ). 

Remarque 

Quand la pulsation d’excitation co est exactement égale à la pulsation propre du 
système, la solution utilisée précédemment n'est plus valide. Le choix de la fonction 
AgCos&V comme solution particulière n'est pas possible car elle est aussi une 
solution de l'équation homogène. La théorie des équations différentielles propose une 
solution particulière de la forme 

x p {t) = t A 0 sin cot 

{cette solution peut aussi être obtenue en faisant tendre co vers C0 Q : vers TD ) 

En reportant cette expression dans l’équation x p + co 2 x p = / cos cot avec 
x p = A 0 sin cot + tcoA 0 cos cot 
x p = coA q cos cot + co A q cos cot -tco 2 A 0 sin cot 
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on obtient 



2 co A 0 coscot = f coscot 

et l'expression suivante pour la solution particulière 

/ 

xit) = t sin eût 

p 2 co 



La solution totale et sa dérivée sont de la forme 



/ 

x{t) = A, shuot + A, cos cot + — tsincot 

2 co 



x{t) = A, cocoscot - A-, cosincot + — sinftV + 

2co 



f tco 

- cos cot 

2(0 



Les conditions initiales 



x 0 = x{0) = A 2 



v, 



0 



= x(0) = A x co 

permettent de déterminer complètement la solution 

V n . / 

sin cot + x n cos cot + 
co 



x{t) = —si 



sin co t + x 0 cos cot + - — t sin cot 
2 co 



La figure 2.2 illustre le phénomène de résonance pour co = CO 0 = y/ k/m quand 
x 0 = 0 et v 0 = 0 . Si le déplacement initial et/ou la vitesse initiale ne sont pas nuis, 
une oscillation libre viendra ce superposer à cette réponse forcée. L'accroissement 
constant de l'amplitude avec le temps montre le caractère artificiel de cette 
modélisation qui ne prend pas en compte l'amortissement. 




Figure 2.2 - Réponse du système non- amorti à une excitation harmonique dont la 
pulsation correspond sa pulsation naturelle ( phénomène de résonance). 



1.2- Notation complexe 

Résultant d'une double dérivation, l'accélération se trouve pour un mouvement harmonique 
en opposition de phase avec le déplacement et s'écrira x(t) = -co 2 x(t). Pour un système 
non amorti, la solution particulière du déplacement sera donc toujours en phase avec la 
force d'excitation en dessous de la pulsation naturelle du système et en opposition de phase 
au dessus. Dans le cas des systèmes amortis, l'équation du mouvement devient 

mx (t) + cx(t ) + kx(t) = Fc os cot 
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Le terme lié à l'amortissement introduit une composante en quadrature (proportionnelle à 
la vitesse) qui se traduit par un retard de phase du déplacement par rapport à la force 
d'excitation. La réponse d'un système amorti (solution particulière) est ainsi une fonction 
harmonique à la même pulsation que l'excitation mais avec une phase différente : 

x p {t) = A 0 cos (cot - </>) 

La même méthode de recherche des solutions que pour le cas non amorti peut être 
employée 1 . Cependant, la méthode complexe permet d'alléger considérablement les 
calculs. Si on écrit pour le second membre de l'équation du mouvement 

F(t)=Fc08a)t = Re{Fe j, °'} 

la solution particulière peut se mettre sous la forme 

x„(î)= Re{Xe J “} 

avec l'amplitude complexe X = A 0 e~ s$ qui contient l'information de phase (le signe 

moins introduit devant la phase met en évidence le retard de phase). Vitesse et accélération 
peuvent donc s’écrire 

x p (t)=Re{jû)X e jC0 ‘} et x p (t)= Re{- co 2 X e Ja '} 

En omettant Re{ }, l'équation du mouvement devient 

(- co 2 m + j co + k)x e jû>t = F e jû> ' 
ce qui permet d'obtenir facilement la solution particulière complexe 

Z7 

Y = A e~ i<p = - 

o e 2 ■ 1 

-co m+ jco + k 

Pour écrire x p (t) on calcule le module et la phase de X 



A 0 = \X\ = ^Re 2 {x}+ Im 2 {x} 



et 



(f) = - arctan 



Im{X } 

Rc{x} 



Par exemple, en considérant la forme ( u complexe) 

u = a + jb et |m| = = ^(a + jb)(a - jb) = V a 2 + b 2 

b 

ç = arctan — 
a 

Cette démarche est illustrée par la section suivante. 



1 e~» 

— = — — j— ou 

u \u\ 

et 



1 par exemple dans D.J. Inman, Engineering vibrations, Prentice-Hall, 1996. 



32 




Vibrations - Acoustique 1 



I I - Vibrations forcées à 1 ddl 



Remarque 

Une variante consiste à utiliser des variables complexes F = F e ia ' et x p = X e ja ' 
(donc tel que F(t)= Rc{f} et x {t)= Re{x /; }). Vitesse et accélération s'écrivent 
x p = jcox p et x = -tir x p et l'équation différentielle 

x + 2 Çco Q x + co g x = — 
m 

1.3 - Excitation harmonique des systèmes amortis 

Sous la forme normalisée par rapport à la masse, l'équation du mouvement en excitation 
forcée s'écrit 

Fit) 

x{t) + 2 Çco Q x{t) + col x (0 = 

m 

avec CO 0 = yjk/m , Ç = c/(2mco 0 ) . En notant comme précédemment la force et la 
solution particulière sous la forme 

F{t) = R e{Fe ja ' } et (t) =Re{v e iat } 

on obtient l’équation différentielle 

(- co 2 + j2Çct) 0 û) + col ) x = — 

m 

et à partir de cette dernière relation, l’amplitude complexe de la solution particulière 



avec 




F jm 

co 2 + j2Çco 0 co 




F jm 

V(®o - co 1 ) 2 + (2 Çco^cof 



A a 



2Çeo 0 co 

et cp = arctan — 

co- - co~ 



La solution particulière x p {t) 




<f) 



est donc de la forme 



x 



.« = 



F !m 



C 



^(col -C0 2 f +(2Çco <] co) : 



: COS 



coi - arctan 



v 



2 Çco 0 co 

2 2 

COq - CD 



La solution générale est la somme de cette solution particulière et de la solution générale 
du système homogène. 

Pour le cas sous-amorti (0 < <r< i) 

x{t) = A e ~^ co ° t sin {co a t + 6 ) + A 0 cos {cot - cj)) 
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où A 0 et (j) sont les coefficients qui viennent d’être déterminés et A et 0 sont obtenus à 
l’aide des conditions initiales (on distingue bien la pulsation naturelle CO 0 et la pseudo- 
pulsation co a ). 

Pour les valeurs importantes de t, le premier terme disparaît et la solution totale 
correspond à la solution particulière : c’est la réponse stationnaire (le premier terme 
correspond à la réponse transitoire). 




Figure 2.3 - Solution générale du système sous-amorti à une excitation harmonique : 
réponse transitoire (ou libre : solution générale du système homogène) et réponse 
stationnaire (solution particulière) 



A et 6 dépendent de x () et v 0 comme pour le système libre mais aussi de F : les 
coefficients sont donc différents du système libre. 

Dans un grand nombre de cas on s’intéresse à la réponse stationnaire du système, sauf dans 
le cas où le système est soumis à une excitation par choc. 

En factorisant par (û\ = k/m , le module X s'exprime sous la forme 



X 



G)- F/k 



V(®o ~^ 2 ) 2 +(2 Ç(O 0 co ) 2 



■ *stat D > 



avec 



le déplacement statique 
l'amplification dynamique D 

réduite r = œ/ (O 0 . 

La phase s'exprime alors par (j) = arctan 



-Xtat = F / k 



(l-rf + (2(rf 
2 (r 



avec la fréquence 



1 - r" 
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fréquence réduite fréquence réduite 



Figure 2.4 - Réponse du système sous-amorti à une excitation harmonique. 



Pour r = 1 [co = (O 0 ) , c’est la résonance et la réponse dépend de P amortissement 



D = 



2<r 



et 



</> = 



n 



mais quand r « 1 alors D ~ 1 et quand r » 1 alors D ~ — . En utilisant des échelles 

r~ 

logarithmiques pour l’amplitude et la fréquence, on obtient la représentation de la figure 



2 . 5 . 




Figure 2.5 - Réponse forcée d'un système à un degré de liberté en représentation logarithmique. 
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Il est important de noter que par principe la résonance apparaît par définition quand 
CO = OJ {) (fréquence naturelle du système non-amorti), ce qui correspond à une rotation de 

phase de 90°. Toutefois, pour un système amorti û) 0 ne correspond pas exactement à la 
fréquence où la réponse en régime stationnaire est maximale. Ainsi 

Ivl F / k a 

a = — = avec r = — , û) Q = 

V(l-r 2 ) 2 +(2 (rf «» 




La valeur maximale de X apparaît à la valeur r = r max qui correspond à l’annulation de sa 
première dérivée 




= 0 



1 

2 



(l — r ) + (2£>)' [2(1 - r 2 \- 2r) + 8^ 2 r]= 0 



-4r(l-r 2 -2^ 2 ) = 0 



donc r max = yj\-2Ç 2 < 1 

Cette relation n'est valable que pour Ç<\/y[2 . Pour des valeurs supérieures de Ç , la 

courbe ne présente pas de maximums car -J 1-2Ç 2 devient un nombre imaginaire. La 
valeur du maximum est 

|X|| F/k 

IU “ tç^\-C 

La fonction de transfert complexe H(jco) ou H(ûj) se définit comme 



H(co) = j r 

H(co) = ^ 
F 

H{a>) = 4 
F 



1 

m {col - or + j2Ço) 0 a>) 

1 

(k - m(0 2 )+ jcCO 

1 

k (l - r 2 + j2Ç r ) 



avec 




col 



k . 

— et r — col O) 0 . 
m 



La fonction de transfert H{co) peut se décomposer en une amplitude et une phase ou en 
une partie réelle et imaginaire comme le montre la Figure 2.6. 

Si la fonction H{co) est représentée dans un espace à trois dimensions (Figure 2.7a), les 
Figures 2.6c et 2.6d correspondent sa projection sur deux plans perpendiculaires. La 
projection de cette courbe sur le troisième plan (plan de Nyquist) est reportée sur la Figure 
2.7b. Au voisinage de co 0 , le lieu de Nyquist est assimilable à un cercle 
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fréquence réduite r 



fréquence réduite r 



Figure 2.6 - Différentes représentations de la fonction de transfert H{Cü) : amplitude 
(a) et phase (b) ou partie réelle (c) et partie imaginaire (d). 



Paramètres du cercle dans le plan de Nyquist : 



Centre : 



f i ( , i+f ! '| i f-i+fV 2 ' |] 

1h~M < Jj 



Diamètre : — 

k 




partie réelle de H(co ) 




(a) (b) 



Figure 2.7 - (a) Fonction de transfert dans un espace de représentation à 3 
dimensions, (b) projection de la fonction de transfert sur le plan de Nyquist 
( Re{//}. Im {//} ) 
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Les pulsations aux lieux des tangentes verticales sont ( 0 ] et co 2 auxquelles correspondent 
les amplitudes \kf{co t )| et \h(co 2 )| 

= « m ax (1 -C) et \H((O x )\ = \H{û) max )|/V2 

«2 = ®,nax (1 + C) et \ h { co 2 )\ = | H (ft» max y V2 

°ù « max = coJl-2Ç 2 est la pulsation correspondant à l'amplitude maximale. 

Ces caractéristiques permettent de déterminer le taux d’amortissement à partir du diamètre 
du cercle ou encore par la méthode de la ‘bande passante à - 3dB’ présentée en fin de 
chapitre (Application 5). 



1.4- Excitation par la base 



Souvent une structure est excitée par l’intermédiaire des plots de suspension (machine 

excitée par les supports, automobile excitée par la route par l’intermédiaire des 

suspensions) , , 

m x[t) 



A x' 




(0 



A 



y(t) 




support vibrant k ( x ?) c (* 

Figure 2.8 - Excitation d'une masse suspendue par la base (support vibrant ) 

En sommant les forces qui s’exercent sur la masse m , on obtient l’équation 



m x(t) + c [x(t) - y(t )] + k [x(/) - y(t )] = 0 

qui peut se mettre sous la forme d'une équation différentielle avec second membre 

m x(t) + c x(t ) + k x(t) = c ÿ(t) + k y(t) 
x(t) + 2Çco (l x(t) + col x ( ! ) = 2^ 0 >’(/) + col y{t) 

On suppose que le support possède un déplacement harmonique 

y(t)= Y cos cot = Re{y e jnn }, Y amplitude réelle 

et on cherche une solution particulière de la forme 

x p (t) = A 0 cos (cot -cf>)= Re{x e JÛ> ' } avec X = A q e J<p amplitude complexe 

En reportant cette expression dans l’équation différentielle précédente 
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(- co 2 + j2Ço) 0 a> + col) X = (j^co 0 co + col) Y 
il est possible d'obtenir l'expression de l'amplitude complexe du déplacement de la masse 

co 0 (û) 0 + j2Çco) Y 



X=A 0 e~ 



col -co 2 + j2Çco 0 co 



a 



Pour calculer A 0 e t </>, on pose X = A 0 e ^ = — (ex ,j8 complexes) 



P 



A> = * = 



f a ^ 



a 

Jj 



xP J 



\a\ 

W\ 



= 7-Z7 = COqY 



^col+{2Çco) 2 



A 0 e ~ JV = — 



a \à\e" a 



X e J 0/1 donc <!> = </>« -<f> a =arctan 



iJ(col -co 2 ) 2 +(2Çco 0 co) 2 

2 ÇcüqCü 



p \p\ 



/? Ya 



-> 2 

CO- - co~ 



- arctan 



2 Çco 
co n 



soit 



X 



,{t) = û) 0 Y 



col + (2 Ccof 



(col - co 2 )' + (2 Cco 0 co) 2 



cos 



{cot - (j)) 



On a vu que la solution particulière pouvait représenter seule la solution stationnaire. En 
utilisant la fréquence réduite r = col co 0 , on appelle transmissibilité en déplacement le 
rapport des amplitudes des déplacements 



PL i 2 

M |_(l-rL + (2<TrL 



Ce rapport décrit comment le mouvement se transmet de la base vers la masse, en fonction 
de la fréquence réduite r (voir Figure 2.9). La transmissibilité est inférieure à 1 quand 
r > ^2 . A la résonance elle dépend de Ç . 




Figure 2.9 - Rapport de transmissibilité en déplacement en fonction de la fréquence 
réduite et pour différentes valeurs du rapport d’ amortissement. 
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1.5- Excitation par déséquilibre dynamique en rotation 

Les machines tournantes constituent des sources de vibrations très courantes. De petites 
irrégularités dans la distribution des masses des parties en rotation causent des niveaux 
vibratoires importants. On schématise une machine de masse m comportant une masse m 0 
en rotation à une distance / de son centre. Un guidage sans friction autorise seulement un 
mouvement dans la direction x (voir Figure 2.10). En supposant la vitesse de rotation C0 R 
constante 

x R (t) = lsvneo R t 



machine 




x{t) 

A 

— > y(t) 



Figure 2.10 - Excitation d’une machine suspendue par une masse en rotation 

La force de réaction générée par la rotation de la masse a une composante dans la direction 
x qui est proportionnelle à m 0 et à l'accélération x R . Cette force agit sur la masse m de 
la machine (les forces dans la direction y ne sont pas considérées). L’équilibre des forces 
par rapport au référentiel du support s’écrit (m 0 fait partie de la masse m de la machine) 

(/ m - m 0 ) x(t) + m 0 [x(t) + x R (t)] = —kx(t) - cJc(t ) 

ce qui conduit à l’équation 

mx(t) + cx(t ) + kx(t) = -m 0 x R (t) . 



Le terme au second membre est la force F(t) = —m 0 x R ( t ) qui excite le système. Ce dernier 
peut s'écrire en normalisant l'équation par rapport à la masse 



x(i) + 2 Çû) Q x(t) + col = 

m 



avec la force qui s'exprime par 



F(l) = 



cl 



(?) = -/ m 0 — — (sin co R t) = l m 0 Cù\ sin û) R t 



-m 0 x R yij--nn 0 

dt 



La méthode d’obtention de la solution particulière est semblable à celle des problèmes 
précédents en utilisant la notation complexe, pour la force 



Fit) = 1 m Q coj, s i n 



Re{-y7 



m 0 col e ^ C ° Rt 



}=Re{Fe Mf } 
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donc l'amplitude complexe est F = -jlm () 03 R et pour l'expression de la solution 
particulière recherchée 



, ( î ) = Aq cos [co R t -</>)= Re{x e j( ° Rt } avec X = A 0 






ce qui conduit à 



ou, avec r = û3 R jo3 () 



X 



- col + 



jm 0 lO)l 

j2Çû) (} û)r ) 



I I /WnZ 

A„ = x - » 



m 



(l-Cf + (2 Çrf 



et 



</> = - 



Imix} 1-r 2 2 Çr 71 

arctan — - r — T = - arctan = arctan - 



Re{X } 



2Çf 



1-r” 2 



Finalement la solution particulière s'écrit 



x p (?) = Aq cos (œ R t -</>) = A q cos 



c 



2 Çr K 



co R t- arctan + — 



V 



f 



1-r” 2 J 



= -Aq sin 



co R t - arctan 



V 



2 Çr 
1 - r 2 ) 



L’amplitude de déplacement en régime stationnaire est une fonction de la vitesse de 
rotation (Figure 2.11). Quand r»l, la valeur normalisée du déplacement tend vers 1 
indépendamment de Ç . Pour les machines tournantes, la suspension sera calculée pour que 
03 () soit en dehors de la gamme de fonctionnement. 




fréquence réduite r 



Figure 2.11 - Réponse du système à 1 ddl à une excitation par une masse tournante. 
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[APPLICATION 1] Accéléromètre 



Une importante application de l’analyse des vibrations forcées et de l’excitation par la base 
se rencontre dans la réalisation de transducteurs pour la mesure des vibrations. C’est le 
cas de l’ accéléromètre. 




m x(t) 



m 



t t 

k{x-y) c{x-ÿ) 



Figure 2.12 - Schéma de principe d’un accéléromètre. 



L’équilibre des forces appliquées à la masse m conduit à la relation 

mx = -c (x - y) — k(x— y) 



La base est solidaire de la structure dont le déplacement est y(t) = Y cos eût = Re{y e Jeo1 }, 
( Y amplitude réelle). Le mouvement de la masse par rapport à la base dont dépend le 
signal électrique produit par le cristal piézo-électrique est 

z{t) = x{t)-y{t) 



L’équation différentielle peut être écrite en terme de déplacement relatif 

m(z + ÿ) + ci + kz = 0 

mz + cz + kz = -tnÿ = mû) 2 Y cos eût 



La solution a la même forme que précédemment 



:« = 



co Y 



e 



^(co 2 -eo 2 f +(2Ça) 0 eo) : 



r cos 



eût - arctan 



v 



2Çû) 0 co 

2 2 

COq - CO 



Le rapport Z /Y des amplitudes correspond à la réponse en fréquence du 



(l - r 2 y + (2 C rf 



avec r = ■ 



m 

Cûn 



\ 

J 

transducteur 



alors que sa réponse en phase est 



(f) = arctan 



2Ç r 
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Fonctionnement en séismographe 

Pour des valeurs importantes de r (r>3), l’amplitude du rapport |Z/F| tend vers 1 : le 
déplacement relatif et le déplacement de la base ont la même amplitude. Il est donc 
nécessaire que la fréquence de résonance C0 0 soit petite car la gamme de mesure débute à 

ttfj ~ 3ty 0 . 

Fonctionnement en accéléromètre 
La solution peut aussi s’écrire sous la forme 



r 2 Y _û) 2 Y 1 

V( i_r 2 ) 2 +( 2 C r f 0)0 V( i_f 2 ) 2 +( 2 C r f 



On remarque que (O 2 Y = Ÿ est l'amplitude de l'accélération et quand 

V(l _r ") 2 + (2Çt'Y — >1, l'amplitude |z| est donc proportionnelle à l’accélération de la 
base Ÿ : 




Ÿ 





Figure 2.14 - Courbes de réponse et sensibilité des accéléromètres en fonction de 
leur masse. 
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La sensibilité est donc inversement proportionnelle à û)q = k/m . La masse m doit être 

faible (i) pour ne pas charger la structure, (ii) pour que la fréquence de résonance soit 
suffisamment élevée (typiquement, plusieurs dizaines de kHz), mais la sensibilité diminue. 

La sensibilité des accéléromètres piézo-électriques est exprimée en pC/ms' 2 ou en pC/g. 
Un conditionneur spécifique (amplificateur de charge) est nécessaire. Certains 
accéléromètres possèdent une électronique intégrée qui réalise le conditionnement du 
signal et leur sensibilité est exprimée en mV/ms 2 ou en mV/g. Un conditionneur au 
standard ICP est alors nécessaire pour alimenter l'électronique intégrée par un courant 
constant (2 mA) par l'intermédiaire du câble signal. 





Figure 2.15 - Alimentation et conditionnement des accéléromètres (en haut ) et 
principaux types d’ accéléromètres : compression (gauche) et cisaillement ( droite ) 

M : masse, P : éléments piézo-électriques, R et S : systèmes de précontrainte 



Sensibilité transversale des accéléromètres 

Les accéléromètres peuvent présenter une sensibilité maximale dont la direction ne 
correspond pas exactement à leur axe. La composante perpendiculaire est une quantité 
parasite qui est quantifiée par sa sensibilité transversale 





Figure 2.16 - Sensibilité transversale 
d’un accéléromètre 
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[APPLICATION 2] Isolation vibratoire : Techniques de base 

L'isolation vibratoire en basse fréquence concerne deux types de problèmes différents : 

° Isolation vibratoire d'une machine : limiter la transmission à l'assise des efforts dus au 
fonctionnement d'une machine, 

D Isolation vibratoire d'un équipement : limiter la transmission du mouvement de l'assise 
à l'équipement. 

A2.1 - Transmissibilité 

Isolation vibratoire d'un équipement : Transmissibilité en déplacement 




Figure 2.17 - Equipement isolé d'une assise vibrante illustrant la transmissibilité en 
déplacement 



La masse est excitée par les vibrations de la base et la transmissibilité en déplacement a 
été précédemment définie dans la section 1 .4 par 



T - 


X 


V 


1 + (2<T rf 


' F 


J 


(i- 


r 2 ) 


2 + (2Çrf 



Isolation vibratoire d'une machine : Transmissibilité de la force 



\ F(,) 




Figure 2.18 -Machine vibrante isolée pour ne pas transmettre de vibrations illustrant 
la transmissibilité de la force. 



45 








I I - Vibrations forcées à 1 ddl 



Vibrations - Acoustique 1 



Sous une excitation stationnaire F(t) = Re{F e jox } le déplacement du système sous-amorti 
est x(t ) = Re{x e J0 * } (ici X = \ e~ J</> est une amplitude complexe), avec 

ii \ F \ 1 

k V(l-r) ; +(2f,f 

La force transmise à l’assise est F t (t) = F k (t)+ F c (t) = Reji 7 , e JÜ *j , soit en notation 
complexe F t - F k + F c 

- la force transmise par la raideur est F k (t) = k x(t) soit F k = kX 

-la force transmise par l’amortissement est F c (t) = cx(t ) soit F c = jcoeX = jû)(2m,û) Q Ç)X 
d'où 

F t = X{k + j2mûW 0 C ) = kX{\ + j2Çr ) et |F t | = k\X\*Jl + {2Çrf 

La transmissibilité de la force est définie par 



T = \ 


If 


L 1 l + (2^r) 2 


F 


\(l-r 2 } + {2Çrf 



Les transmissibilités de la force et du déplacement correspondant à deux problèmes 
différents ont cependant la même expression. La transmissibilité peut également être 
exprimée en dB : 



L t = 201og 10 T [dB] 



La transmissibilité est représentée en linéaire sur la Figure 2.19a et en dB avec une échelle 
logarithmique des fréquences sur la Figure 2.19b. 





Figure 2.19 - Transmissibilité en fonction de la fréquence réduite et pour différentes valeurs du 
rapport d’amortissement, a) représentation linéaire (gauche), b) représentation en dB (droite). 
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Remarques 

1) Pour différentes valeurs du taux d'amortissement Ç les courbes passent toutes par le 
point T = 1 et r = V2 . 

2) Pour que la force transmise à l'assise soit inférieure à la force d'excitation ou que le 
déplacement de la masse soit inférieur à celui de la base il faut que la fréquence réduite 
soit 

(0 f rr 

r = — = — > V2 
(Oo fo 

3) Dans la zone d'atténuation (r>V 2), la transmissibilité est pratiquement 
indépendante du taux d'amortissement Ç si celui-ci est inférieur à 10%. Elle peut alors 
s'approcher par la formule 




Quand r » 1 , la transmissibilité exprimée en dB vaut approximativement 
L T ~ -lOlog r , c’est à dire qu'elle décroît de 40 dB/décade (r 2 /rj=10) ou de 
12 dB/octave ( r 2 / r, = 2 ). 



A2.2 - Atténuation 

Pour les opérations de réduction du niveau vibratoire, le taux d'atténuation (souvent 
exprimé en %) est également utilisé 




Attention : Quand l'atténuation (ou gain ) est exprimée en dB, c'est la quantité 
20 log 10 — = —L r ( positive pour r > ) qui est employée. 



La fréquence d'excitation / pour laquelle on obtient une atténuation A se déduit des 
équations précédentes : 

= IÏEÏ 

a 0 fo VA — 1 

Pour une fréquence d'excitation / , l'atténuation A dépend seulement de la fréquence 
propre / 0 de la suspension. Ces résultats sont présentés sous forme d'abaques (Figure 
2 . 20 ). 
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1 2 

10 10 

fréquence d'excitation f 




fréquence d'excitation f 



Figure 2.20 - Abaques d'atténuation vibratoire donnant les relations entre la fréquence d'excitation 
f et la fréquence propre f 0 . a) pour des atténuations A = \ — T , b) pour des affaiblissements 

Att = 201og 10 (l/r) (en dB). 



A2.3 - Détermination d'une suspension élastique 

Un cas fréquent est celui d'une excitation due à un mécanisme en rotation. Une atténuation 
A est recherchée dans une plage de fonctionnement définie par l'intervalle de fréquence 
[fi ’f 2 J • C’est la fréquence d'excitation minimale f l qui représente la contrainte la plus 
importante et qui sert à déterminer la fréquence propre f 0 de la suspension élastique. 

Exemple : une atténuation A = 90%, correspond à T = 1 - A = 0. 1 , soit un gain de 
20 dB. Si la plus petite fréquence à atténuer est f l = 30 Hz, l'abaque donne une 
fréquence propre de la suspension de 9 Hz (qui se calcule aussi par 
/ 0 = / A /(A-l)/(A-2) = 30^(0. 9 - 1)/(0.9 - 2) = 9 Hz. La relation û) 0 = g /x 0 

permet de connaître la déflexion statique : x 0 = g/4;r 2 / 0 2 =9.8l/4/r 2 9 2 =0.003 m. 

Le centre de gravité de l'ensemble mécanique à suspendre permet de déterminer le nombre 
et la position des points de fixation pour que la charge soit également répartie (dans le cas 
où la position des points de fixation est imposée, on peut être conduit à choisir des supports 
de rigidités différentes). 
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Exemple : Pour une charge de 300 daN répartie sur 6 points, il faut chercher dans le 
catalogue du fournisseur des supports ayant une fréquence propre de 9 Hz (ou une 
déflexion statique de 3 mm) sous 50 daN. 

Dans le cas d'une machine tournante, le régime de fonctionnement est atteint après le 
démarrage en passant par la fréquence propre ( r = 1 ). La transmissibilité passe par un 
maximum qui est contrôlé par le taux d'amortissement 

T = — 

max 

Les élastomères utilisés dans les suspensions ont des taux d'amortissement compris entre 
0.02 et 0.1 et jusqu’à 0.2 pour des synthétiques. 

Exemple : Au passage à 9 Hz, le support qui possède un taux d'amortissement de 5% 
(£ = 0.05) assure une transmissibilité 7j mx =10, soit 20 dB. Pour une machine 
suspendue, l'amplitude du déplacement est X rnax ~ {E jk)T mm . Pour un équipement 
suspendu, l'amplitude des vibrations de l'assise est multipliée par 10. 



Remarques 

1) la plupart des supports en élastomère de formes géométriques simples ont un 
comportement "presque" linéaire dans la plage de fonctionnement recommandée. Par 
exemple la raideur est donnée par (voir Ligure 2.21) : 



support en compression : 
support en cisaillement : 



k = 



ES 

~T 

GS 

h 



(Ligure 2.21a) 
(Ligure 2.21b) 



avec 



E module d' Y oung [N/m 2 ] 

S section du support 

/, h épaisseur du support 

G = E/ 2(1 + v) module de cisaillement ( v coefficient de Poisson). 




Figure 2.21 - Supports élastiques travaillant : a) en compression, b) en cisaillement. 
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Figure 2.22 - Différentes formes de supports anti-vibratoires. 



2) dans le cas de support dont la rigidité statique diffère de la rigidité dynamique, la 
relation entre la fréquence propre de la suspension et l'écrasement statique n'est plus 
valable. Par exemple, les supports "Evidgom" de Paulstra ont une caractéristique de 
rigidité représentée sur la Figure 2.23. La rigidité dynamique est déterminée par la 
tangente au point de fonctionnement. 



tangente 



Charge * k 




Figure 2.23 - Rigidité dynamique et rigidité statique d'un support élastique ("Supports élastiques", 

Documentation Paulstra ). 
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2 - DIFFERENTES FORMES D’AMORTISSEMENT 

2.1 - Amortissement par frottement solide (ou amortissement de Coulomb) 

L'amortissement de Coulomb est l'amortissement par frottement solide. Il est difficile de 
donner une forme mathématique correcte pour modéliser un amortissement. Contrairement 
à la masse et à la raideur, l'amortissement ne peut pas être déterminé par des tests statiques. 
L'amortissement visqueux en ex est la forme la plus courante mais l'amortissement par 
frottement solide se rencontre fréquement dans les mécanismes. Il est défini par la relation 

- juN x>0 

F d= F d(x) = < 0 * = ° 

/iN x < 0 

F d force de dissipation, 

N force normale, 

fl coefficient de frottement. 

Figure 2.24 - Frottement solide : notation 

La force F d est toujours opposée à la direction du mouvement (son signe est contraire à 
celui de x ) 



x{t) 




mx + kx= /uN pour x < 0 

jUN = jUtng 



mx + kx = -jUN pour x > 0 



-jU N = -jumg 

Figure 2.25 - Décomposition de l'équation du mouvement pour un frottement solide. 
Dans ce cas particulier, la niasse est posée sur le plan et N = mg . 

C'est le signe de x qui détermine la direction de la force de frottement 

mx + juN sign(x) + kx = 0 

avec la fonction 'signe' 
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sign(w) = 



' 1 
< 0 
-1 



quand u > 0 
quand u = 0 
quand u < 0 



Cette équation ne peut pas se résoudre directement par les méthodes employées 
précédemment car c'est une équation différentielle non-linéaire. Pour obtenir la réponse 
libre du système, on va considérer séparément les tronçons temporels entre les 
changements de direction de x . 

Conditions initiales : x(0) = x 0 et x(0) = 0 

x 0 se trouve à droite du point d’équilibre et la force kx 0 est plus forte que F d 




i (0 ) = 0 



Déplacement vers la gauche : x<0 

La masse se déplace vers la gauche et l'équation valable pour x < 0 a une solution de la 
forme 



avec co Q 
initiales 



t \ > r, ■ MM 

x\t ) = Aj cos co 0 t + B x sin co 0 t -I 

k 

yjk/m . A l et B x sont des constantes qui sont déterminées à l'aide des conditions 



*(o) = A +^- = x 0 

k 

x(0) = =0 



a -OÎL 

k 

B x = 0 



Le mouvement vers la gauche à partir du point de départ x 0 est donc 



t(f) : 



jUN 



cos co 0 t + 



pN 



et il se poursuit jusqu'à l'instant /j où la vitesse x = 0 



x{t l ) = -û) o 



pN 



sin co 0 t 1 =0 



K 



c'est à dire à t x - — . Ensuite, à condition que kx > [A ’N , la masse se déplace vers la 



co. 



droite. 



Déplacement vers la droite : x > 0 

L'équation différentielle correspond au cas où x > 0 et a pour solution 
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, x „ . UN 

c{t ) = A 2 cos co Q t + B 2 sin co Q t 

k 



n 



pour — <t <t 2 . Les conditions initiales pour cette équation sont fournies par l'équation 



t o n 



précédente à t 1 



ce qui conduit à 



f n ^ r 



v®oy 
f n ^ 

j 



jUN 



v 0 



V k j 



COS 71 + 



jUN 2juN 



■X n 



= -co , 



x, 



jUN 



0 



V 



sin 7l = 0 



) 



A 2 x 0 



3 jUN 



La solution pour le second intervalle de temps est 



B 2 = 0 



c(r) = 



x n 



3/J N 



cos û) 0 î ■ 



jUN 



n 2 n 

— < t < — 



OA 



co n 



Arrêt du mouvement 



Cette procédure se répète jusqu'à ce que le mouvement s'arrête. Cet événement se produit 
quand 

x(t) = 0 et |fcr| < /J N 




temps [s] 



Figure 2.26 - Oscillation libre avec frottement solide 
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Remarques : 

a) L'amplitude des oscillations décroît linéairement comme 

2/u Nco t 

x o ; ^ 

7tk 

b) La fréquence d'oscillation du système amorti par frottement solide est la même que 
pour le système non amorti. 

c) La position d'arrêt peut être différente de la position d'équilibre initial. 



2.2 - Energie dissipée par cycle et amortissement équivalent 

La réponse harmonique forcée d'un système avec amortissement de par frottement solide 
est la solution de l'équation 

mx + fi N sign(x) + kx = F cos cot 

Plutôt que de résoudre cette équation directement, on recherche une solution approchée en 
considérant la solution du système à amortissement visqueux qui dissipe la même énergie 
par cycle. C'est une hypothèse raisonnable tant que la force d'excitation est bien plus 
importante que la force de friction (F » jliN ). 



On suppose une réponse stationnaire de la forme 

x(t) = X sin cot 

L'énergie dissipée par cycle A E par un système dont le coefficient d'amortissement 
visqueux vaut c est 



AE V = f , F d dx = 

J cycle 




avec x{t) = coX cos cot 



l7tj(0 271/(0 

A E v = J cco 2 X 2 cos 2 cot dt = J -co 2 X 2 dt = nccoX 

o o ^ 



, 2 1 1 „ , 

(car cos cot = — + — cos2<yt ) 

2 2 

Par ailleurs, l'énergie dissipée par cycle par frottement solide est 



27tjû) 



A E s =^ F d dx= J jU ’N sign(i)— dt=juN j sign(x) x dt 



dt 



Ik/co 



En séparant cette intégrale en 3 zones en fonction du signe de x 
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A E s = juN coX 



7tj2co 2>7 tI2cq 2k jœ 

J cos cot dit- J cos cot dt+ J cos cot dt 



k/2co 



B 



2Kj2(û 



C 



K/2CO 



= 4/u N coX |cos<wr dt = 4juNX 

o 

1 — 1 — 1 

(en considérant — \ 2 cosudu = — [siru/l 2 = — ). 

m JO m L -*0 ni 




Pour qu'un système à amortissement visqueux dissipe la même énergie qu'un système à 
amortissement par frottement solide : AE V = AE S , il faut qu'il ait un coefficient 
d'amortissement visqueux équivalent 

_ 4 ju N 
eq 71(0 X ’ 

ou en terme de rapport d'amortissement équivalent 



c eq _ 2/jN /m 
eq 2 mo 0 7T(O 0 (O X 



c q = 



On constate qu'il dépend de l'amplitude et de la fréquence d'excitation. Ainsi le système à 
frottement visqueux dont le mouvement est décrit par 



P 

x + 2 C m û) n x + col x = — cos co t , 
q m 



dissipera la même énergie par cycle que le système à frottement solide. Sa solution est de 
la forme x{t) = X cos (cot -(/>), X amplitude réelle 



X = 



F/k 



F/k 



l(l-r 2 f+(2C„ rf ^(l-r 1 ) 2 +(4^N/xkXf 
En résolvant cette expression pour X 



X 



(i -sr+4 



x 



f 4 jI/N) 


2 " 






{ TTk J 






l k) 



on obtient pour l'amplitude 



X 



Vl-(4 juN/ttF) 2 

I l-r 2 | 



Par ailleurs, la phase s'écrit 



l = arctan EL* r = arctan 4> f N 

1 - r 2 7rkX{\-r 2 ) 
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En substituant dans cette dernière équation l'expression de X 



±AuN 

ç = arctan 


(le signe ± dépend de 1 - r 2 


iïF-yJl - (AjUN /7rF) 2 





ainsi quand r < 1 => (j) > 0 et quand r > 1 => (j)< 0 . 

Remarque 

Contrairement au système à amortissement visqueux, à la résonance ( r = 1 ), le modèle de 
système à frottement solide présente une amplitude infinie et une discontinuité pour la 
phase. 

Les expressions pour la réponse sont valables seulement dans le cas où AjuN < K F , ce qui 
confirme l'hypothèse de départ. 

L'amplitude du système n'est pas une fonction linéaire de la force d'excitation et sa phase 
dépend de F . Il ne satisfait pas au principe de supeiposition. C'est un système non- 
linéaire. 



2.3 - Amortissement hystérétique (ou structural) 

Nous avons vu qu'il est intéressant d'étudier les mécanismes d'amortissement en examinant 
l'énergie dissipée par cycle définie par l'équation 

À E = f F d dx, F d : force d'amortissement 

J cycle ( 

On définit une capacité d'amortissement spécifique par le rapport de l’énergie perdue par 
cycle sur la valeur crête de l'énergie potentielle AE/U ^ . Le facteur de perte est 
l'expression précédente normalisée par 2 n 

A E 



U uax est l'énergie potentielle au maximum du déplacement. 



Pour un système à amortissement visqueux 

kccoX 2 _ ceo 
71 ~ 27i(\kX 2 )~ k 



A la fréquence de résonance a> = û) 0 = yjk/m , le facteur de perte est le double du rapport 
d'amortissement : 



V = 




2 Ç. 
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Si on considère la force nécessaire pour déplacer la niasse d'un système à amortissement 

visqueux 



F = kx + ex 



avec x = X sin cù t et x = (û X cos Cû t . Avec cos co t = ±^l - sin 2 cot , il est possible d'écrire 

F = kx± c(Oyl X 2 - x 2 

Dans le plan (F,x), cette relation est représentée par l’ellipse de la Figure 2.27a. C'est une 
boucle d'hystérésis dont l'aire correspond à l'énergie perdue par cycle A E v =nc(oX 2 . 
Pour c = 0 , l'ellipse devient la droite F = kx . 




Figure 2.27 - (a) Ellipse représentant la boucle d'hystérésis dans le cas d'un 
amortissement visqueux, (b) Boucle d'hystérésis dans le cas d'un amortissement 
structural. 

En mesurant les contraintes et les déformations de matériaux excités par des efforts 
harmoniques, il apparaît également des boucles d'hystérésis semblables. Elles sont dues à 
des dissipations d'énergie appelées amortissement hystérétique ou amortissement 
structural. L'aire fermée dans la boucle correspond à la perte d’énergie par cycle 

A E h =7rkj3X 2 

qui est proportionnelle à la raideur, au carré de l'amplitude et à /3 la constante 
d'amortissement hystérétique ( fi est indépendant de la fréquence). Si on compare cette 
dernière expression à celle de la perte d'énergie par cycle par amortissement visqueux 

A E v = A E h => nc^toX 2 = nkpX 2 

le coefficient d'amortissement équivalent d'un système à amortissement visqueux qui 
dissiperait ma même énergie par cycle 
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c 



eq 



k/3 

co 



y? étant déterminé expérimentalement à partir de la courbe d'hystérésis. L'équation du 
mouvement d'un tel système est 



..J3k. 

mx H x + kx = b cos cot . 

co 



Si la réponse en régime stationnaire est de la forme x(t)= X cos {cot -(/>), on obtient le 
même type de solution que précédemment 



X = 



F/k 



(l-r 2 )" +(2 Ç r) 



et </> = arctan 



2 Cq r 

1 - r 2 



où Cq = c f . q /2y[km = /3co Q /2oo = (3 /2r , soit 



X = 



F/k 



(1 -rf+f 



et (j) = arctan 



P 



1 -r 2 



Contrairement au cas de l'amortissement visqueux, on observe que l'amplitude maximale 
de Xk/ F à lieu à la fréquence de résonance r = 1. Autre particularité, la réponse n'est 
jamais en phase avec la force d'excitation. 

En utilisant des variables complexes (voir la remarque de la section 1.2) 

mx + + kx = F e jü>t 

co 



avec x=X e JÛ> ' et x= jcoX e jm ' = jcox, d'où 



mx + k(\ + j/3)x = F e ja>t 

L'amortissement hystérétique (ou structural) peut donc être pris en compte sous la forme 
d'une raideur complexe 

k = k{l+ jj3). 
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3 - REPONSE FORCEE A UNE EXCITATION QUELCONQUE 

3.1 - Réponse impulsionnelle 

La réponse impulsionnelle d’un système à un degré de liberté se définie comme la réponse 
à une excitation de la forme F(t) = F 0 £(t) où ô{î) est la fonction de Dirac : 

x{î) = F 0 h(t) 



La réponse impulsionnelle est donc 





o, 


pour t < 0 


h(t) = 


— - — exp (-Çcoj) sin (ût. 


pour t > 0 




mû) n 

l a 





avec la pseudo-pulsation co a = (O 0 y]l - Ç 2 , 

La réponse impulsionnelle est donnée pour une impulsion unitaire Si t ) appliquée à t = 0 . 
Une excitation impulsionnelle 2 appliquée à l'instant t = t l 

F{t) = I 0 S{t-t 1 ) 

aura pour réponse 

x{t)= I Q h(t - t x ) = — — e _ ^ ,_ ^sin (Q, (t - 1,) 
mû) a 

pour t > t x et x(t) = 0 pour t<t x . Cette solution s'entend au sens d'une solution 
particulière, une solution générale viendrait s'ajouter si les conditions initiales n'étaient pas 
nulles (existence d'un mouvement pour t < 0). 

3.2 - Réponse à une excitation non-sinusoïdale 

La réponse à une excitation quelconque Fit) correspond au produit de convolution de 
cette excitation par la réponse impulsionnelle du système à 1 ddl 

+oo 

xit) = Fit) *hit) = J Fir)h(t - t)c1t 
Puisque la réponse impulsionnelle est causale (hit) = 0 pour t < 0) 

t 

xit) = ^ Fir) h{t - t) dr . 

t 

Si l'excitation est-elle aussi causale, on pourra écrire x(t) - J F{t) h{t- r)dr . 

o 



2 L’impulsion J' f’ (t)dt =J" / 0 S(t)dt = I 0 s’exprime en Ns ou kg m s' 1 et l„h(t) a bien la dimension d’un 
déplacement. 
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Figure 2.28 - Calcul de la réponse à une excitation quelconque par convolution par la 
réponse impulsionnelle du système à un ddl. 



Exemple : Réponse de l'oscillateur masse-ressort amorti à une force en échelon (fonction 
de Heaviside) 

f° 

mx + cx+ k x = i _ 

l F o 

Dans le cas d'un système sous-amorti, la réponse s'écrit 

x(t) = [ -^—e~ (c0o(t ~ r) sinû) a (t - z)dr 
l mC °a 

= e'- c ° at f e ?a ° T sin co a ( t - z)dz 
mCO a t 



t<t o 
t>t 0 



L'intégration de cette fonction (table) conduit à 

x{t) = -p - —j = ^ e~^ (t ~ t0l) cos [û) a {t -t Q )-(f>\ pour t > t ( 

r 

avec0 = arctan . - 
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En absence d'amortissement, le système continue à osciller indéfiniment 

J? 

x(t ) = — (l- cos&> 0 ?) 
k 

Quand 0 < Ç < 1 et ?» 0 , on tend vers la réponse stationnaire du système 




La Figure 2.29 représente l’excitation et la réponse dans cette configuration 



F(t ) 2 T T T T T 

1 ' 

o 

-i - 

2 i i i i i 

0 0.5 1 1.5 2 2.5 3 




Figure 2.29 - Excitation (haut) et réponse (bas) d'un système à 1 ddl. 
m = 10 kg, co 0 = 20 rad/s, Ç = 0.05 



3.3 - Réponse à une excitation périodique quelconque 



Soit une excitation qui correspond à une fonction ayant une période T, telle que 
F{t + T) = F{t ) . Cette fonction peut être représentée par une série de Fourier 



F(t) 



+ ^ ( a n cos nco T t + b n sin nœ T t ) , 



n = 1 



L'équation du mouvement est 
et sa solution particulière 



mx + c x + k x = F{t) 



co T 



2n 

T 



X p {t) = X 0 (?) + X x n k COS {nû) T t -(/)„)+ b n sin [nû) T t - <f> n )] 

n = 1 
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avec 



X 



n 




+ (2 Ço 0 nû) T ) 2 




2 Çco 0 nco T 

et è n = arctan 7 -7 

col-{ncü T \ 



La solution particulière x 0 {t) satisfait l'équation 



mx 0 + cx 0 + k x 0 



-> 



X n 



« 0 . 

2k 



Pour obtenir la solution complète, on ajoute la solution de l'équation homogène (la réponse 
libre). Pour le cas où 0 < Ç < 1 , 

x(t) = A sin( co a t — 9) + x p (t) 



où A et 0 sont déterminés par les conditions initiales : le premier terme représente donc la 
réponse transitoire et x p {t) la réponse stationnaire. 



3.4 - Utilisation de la transformation de Laplace 

La transformée de Laplace est un outil particulièrement utile pour calculer la réponse des 
systèmes à un grand nombre d'excitations. C'est aussi une méthode générale pour résoudre 
les équations différentielles. La transformée de Laplace d'une fonction intégrable f{t) 
telle que f{t) = 0 pour t < 0 

£{f(t)}= F(s) = \ f(t)e~ st dt 
0 

5 est une variable complexe. F{s) est l'image de la fonction f{t) dans le plan 
complexe s . La transformée de Laplace possède des propriétés intéressantes pour la 
dérivation 



£{x(t)}= sX (5) - x( 0 ) 
.2?{x(?)}= s 2 X(s)~ .s'x(o) - x(o) 



où v( 0 ) et i( 0 ) sont les valeurs initiales. Il est donc possible de résoudre une équation 
différentielle en tenant compte de ses conditions initiales. Par exemple, le système à 1 ddl 
non-amorti décrit par 

x(t) + col x(f) = 0 x(0) = x 0 et x(0) = v 0 



La transformée de Laplace de l’équation différentielle conduit à 

s 2 X(.sj - sx 0 - v 0 + col = 0 



d'où 



x( s ) 



2 , ,,2 

^ + co Q 



Pour revenir à la description temporelle de la réponse, on emploie la transformée inverse 
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I +°° 

{x[s)}= x(t) = J X[s)e st ds 



Plutôt que de calculer l'intégrale, en pratique on utilise souvent une table ou "dictionnaire 
d'images" : 



x(î) = 



SX n 



■ + ■ 



2 , ,J2 2 . ,J2 

S + CO 0 s + co 0 



— -!> = cosaf et 



2 . 2 

S + Cl 



2 , 2 

S + Cl 



= sin at 



, , Vq 

d'où la solution : x(t) = x n cos COJ + — sincoJ 



ax 



Exemple : Calcul de la réponse impulsionnelle d'un système à 1 ddl sous-amorti 
Pour une impulsion unitaire, l'équation du mouvement est 

mx{t) + cx{t) + k x{t) = ô{t) avec * 0 = 0 et v 0 = 0 

La transformée de Laplace conduit à 

[ms 2 + C5 + k ) X{s) = 1 

ou en normalisant par rapport à la masse 



d'où la solution 



(s 2 + 2Çû) () s + col ) X{s) = — 

m 



, \ V m I \ 

X(s) = 2 —r y = H(s) 

s - + 2 Çco Q s + cOq 



Un dictionnaire d'images (J. Hlaclik, La transformée de Laplace à plusieurs variables, 
Masson, 1969) donne 



1 



p + 2 abp + b~ 



donc 



<-> 



c(f) 



b 2 (\-a-)] 



2 e~ ab ‘ sinl 






■a 2 t 



pour a 2 < 1 



l/m 






e c ° v sinl ù), 



J o 






qui est bien la réponse impulsionnelle du système 

x[t)= £ a {h[s)\= h[t) 



La transformée de Fourier est aussi largement employée 

+oo 

X(co) = \x(t)e- jm dt 



i T 

x{t) = — X{co) e jat cleo 
2 71 J 



63 




I I - Vibrations forcées à 1 ddl 



Vibrations - Acoustique 1 



La transformée de Fourier peut s'obtenir à partir de la transformée de Laplace avec s = jco 
et en changeant les bornes d'intégration 

+oo 

H{cû) = H{jco) = \h{t)e~ jm dt 



3.5 - Réponses à des excitations aléatoires 

Pour les excitations déterministes les fonctions du temps sont connues. Avec des 
excitations aléatoires, la réponse du système est seulement connue sous forme statistique. 
Aussi la réponse x{t) ne pourra être obtenue que sous une forme identique. 

La densité de puissance spectrale (autospectre) est une mesure de la répartition de 
puissance dans le domaine des fréquences d'un signal aléatoire 

S Jco) = lim {X\co, T)X (co, T)} 



h+T 

où x{co,t)= J x{t) e~ jû)l dt 

'i 

est la transformée de Fourier d'un tronçon de signal de durée finie T (pour assurer la 
convergence). La réponse d'un système de réponse impulsionnelle h{t) excité par une 
force aléatoire F{t) est représentée par le produit de convolution 

x{t) = F{t) *h{t ) 



dont la transformée de Fourier est 

X(co, T) = F{F(t ) * h{t))= H(co)F(co, T) 

Fa densité de puissance spectrale de la réponse x{t) s'écrit alors 

S Jco) = H{co)H\co) lim jJf{cd, T)F(a>, T)} 

Sjco)=\H{cofS F Jco) 



Pour déterminer la fonction de transfert H(co) d'un système avec un signal aléatoire (pot 
vibrant excité par un bruit blanc par exemple), on utilise la densité de puissance 
interspectrale (interspectre) entre le signal de la réponse et celui de la force d'excitation 

S Jco) = lim yE{F(co. T)X(ox T)}. 
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En remplaçant X(co,T) par H(a>)F(co,T) on peut écrire 



H(œ) = 



S Fx {a>) 

S ff {oj ) 



Remarque : La fonction de transfert obtenue par cette procédure est connue comme étant 
de type Hi . 



[APPLICATION 3] Analyseur de spectre 

Cet analyseur est très utilisé en vibration et en acoustique. Il permet par exemple de faire 
une estimation de la densité de puissance spectrale ou interspectrale (plus brièvement, 
interspectre ou autospectre). Il est généralement basé sur la transformée de Fourier rapide 
(FFT) d'un signal échantillonné et numérisé. 

Par exemple, l'interspectre entre deux signaux x(t) et y(t) s'exprime mathématiquement 
par 

Sjf)=ümjE{x*(f,T)Y{f,T)} 

Comme cet interspectre représente une distribution de puissance sur une échelle de 
fréquence allant de [- oo. + oo] , il est nécessaire d'en tenir compte pour l'intégrer sur une 
bande [/,,/,] : 

~fl fl 

p= lsjf)df + jsjf)df 

-fl fl 



Pour des signaux x{t) et y(t) réels, l'interspectre présente des symétries (paire pour la 
partie réelle et impaire pour la partie imaginaire) qui conduisent souvent à utiliser une 
définition (interspectre unilatère) qui permet de ne considérer que les fréquences positives 

Gjf) = ]imjE{x-(f,T)Y(f,T)} pour /> 0 

fl 

tel que P=^G xy {f)df puisque G xy {f) = 2S xy {f ) pour /> 0. 

fi 



Remarque : les relations de symétrie pour les interspectres de signaux réels dont il a été 
question plus haut conduisent aux définitions : 

Sjf)=G(f)/ 2 et S (- / ) = G* (/ ) / 2 



En pratique, un tronçon du signal x(t) de durée T est échantillonné selon le critère de 
Shannon (At = 1/ F e , F e : fréquence d'échantillonnage), puis numérisé pour produire un 
vecteur de M valeurs ( T - M At , résolution : Af = 1/ T = F e / M ) : 

x m = x (t m \ m = {\ ■■■ M) 
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La FFT appliquée à ce vecteur conduit à Ml 2 valeurs pour la partie réelle et M 12 
valeurs pour la partie imaginaire : 

i JV . 27tk{m - 1 ) 

X(/»)=FFT{x„}=— 

Dans le cas de signaux aléatoires (bruits par exemple), l'interspectre est estimé avec un 
nombre limité de tronçons de signaux et l'espérance mathématique E est remplacée par une 
moyenne sur N tronçons (N vecteurs de longueur M) 

L'estimation de l'interspectre s'écrit alors 

â»(/*)=|- î x:(/ k )?.(/>) 

* ’ n = 1 

La Figure 2.30 représente le schéma fonctionnel d'un analyseur de spectre pour obtenir un 
interspectre selon l'expression précédente. 




Figure 2.30 — Schéma fonctionnel d'un analyseur de spectre 



Remarques 

1) un filtre anti-repliement permet le respect de la condition de Shannon avant 
échantillonnage du signal (échantillonneur-bloqueur et conversion analogique- 
numérique), 

2) une fenêtre (hanning) est souvent utilisée pour multiplier le signal numérique pour 
réduire l'effet de troncature (création de lobes secondaires), 

3) l'autospectre G xx {f k ) s'obtient en remplaçant y{t) par x(t), 
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4) la formulation de G (f k ) associée aux définitions des transformées de Fourier 
discrètes X n (f k ) et Y n (f k ) conduit à un interspectre dont la puissance de chaque 
raie spectrale vaut 

G,(A)=jG,(/)d/ 

A f 

Pour obtenir une estimation de la densité spectrale à la fréquence f k qui soit 
indépendante de la résolution À/, il faut diviser G {f k ) par Af (s'assurer que 

cette opération est réalisée par l'analyseur en calculant le spectre à des résolutions 
différentes, en changeant M ). 



[APPLICATION 4] Mesure des fonctions de transfert 



On trouvera différentes formes de fonction de transfert en fonction de la réponse 
(déplacement ou accélération) 



Accélérance 



A(co) 



a((ü) 

JW 



avec l’accélération a{co) 



Admittance Y ( û )) = 



x(co) 

JW 



A(c y) 



œ 



avec le déplacement x{û)) = ■ 



a(o>) 

..2 

Cû 



La fonction de transfert mesurée (FRF : Fonction de Réponse en Fréquence) dépend du 
type de capteur utilisé pour mesurer la réponse de la structure. Nous désignons par H Fa (ûj) 
la fonction de transfert mesurée entre un signal de réponse délivré par un accéléromètre et 
un signal d'excitation délivré par un capteur de force. 




Figure 2.31 - Détermination de la fonction de transfert entre deux points d'une structure 
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La Figure 2.31 illustre la détermination d'une fonction de transfert entre deux points d'une 
structure. Deux types d'excitations sont employés : un pot vibrant ou un marteau d'impact. 
Dans les deux cas, le signal d'entrée est délivré par un capteur de force (Fi ou F 2 ) et le 
signal de sortie provient de l'accéléromètre (a). 

a) pot vibrant 

C’est un moteur électrodynamique piloté par un signal aléatoire (bruit blanc) ou 
déterministe large bande (balayage en fréquence, chiip, ...). La force est mesurée 
par un capteur de force placé directement sur la structure. La liaison entre le capteur 
et le pot vibrant est parfois assurée par une tige ( stinger : corde à piano ou nylon) 
dont le comportement en flambage va supprimer l’influence des défauts 
d’alignement du pot qui auraient tendance à exciter la structure par un moment. 



Threaded Insert 



Rubber Dirt 
Shieid 



Tangential 

Flexture 



Magnetic 

Structure 



Rubber Feet 




Table 



Drive Coll 



Terminal 

Radial 

Flextures 



Permanent 

Magnet 



structure 




Figure 2.32 - Pot vibrant et son montage sur une structure avec un capteur de 
force. 



Ce mode d’excitation ajoute localement une masse qui devra être corrigée pour les 
structures légères. La fonction de transfert est calculée à partir de densités de 
puissance spectrale (interspectre G Fa (û)) et autospectre G FF (co) ) 

G ff 

(estimateur de type H1 : le bruit en sortie tend à être éliminé). 
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b) marteau d’impact 

Le capteur de force est placé dans la tête du marteau. C’est une excitation large 
bande dont la fréquence maximale utilisable f c dépend de la raideur de l’embout. 



Force, N 






Figure 2.33 - Marteau d'impact et représentation schématique de l'enveloppe 
spectrale de l'excitation en fonction de la nature de l'embout utilisé. 



Il n’y a pas de problèmes de masse ajoutée. L’énergie transmise à la structure reste 
faible mais la technique possède l’avantage d’être rapide et facile à mettre en 
œuvre. Elle demande cependant une certaine habileté à l’expérimentateur pour 
éviter les rebonds. 

La fonction de transfert est obtenue à partir du rapport des transformées de Fourier 
des signaux d’excitation et de réponse (c'est à dire du rapport G Fa (f k )/G FF (f k ), 
calculer sur un seul tronçon de signal, N = 1), souvent moyennée sur plusieurs 
réalisations. Le signal impulsionnel de force étant de courte durée, une fenêtre 
rectangulaire est parfois employée pour réduire le bruit (et les rebonds ?). La 
réponse d’un système peu amorti risquant d’être tronqué par la fenêtre temporelle, 
une fonction exponentielle doit être utilisée. Cette procédure étant équivalente à 
l’ajout d’un amortissement supplémentaire, une correction devra être apportée par 
le logiciel dans le cas de son utilisation pour la détermination de l'amortissement. 
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[APPLICATION 5] Mesure de l'amortissement 2 

Les propriétés mises en évidence dans les sections 1.2 et 1.3 ont conduit à deux méthodes 
de détermination du facteur d'amortissement (qui viennent s'ajouter à la méthode du 
décrément logarithmique) valables pour les faibles valeurs d'amortissement (^<10%) et 
largement employées en analyse modale expérimentale : la méthode de la bande passante à 
3 dB et la méthode d'identification du cercle. Ces deux méthodes utilisent des fonctions de 

transfert obtenues expérimentalement. C'est l'admittance Y (m) = [ = - 4^1 q U j es t 

F\co) 0)- 

employée, mais on voit comment elle peut se calculer facilement à partir de la fonction de 
transfert accélération/force H Fo (ûj) = A(co). 



Méthode de la bande passante à 3dB 
Le module 

\y(oi\=- 






m 



présente une valeur maximale à la pulsation &> max = to oy j 1 - 2Ç 1 (voir section 1.2). Poul- 
ies faibles amortissements que l'on rencontre dans les structures métalliques (inférieur au 

I T 1 

1%), ûJ nmx = co, et 7 max = m VU) - «L ) 2 + (2^ 0 «max Y =l/(2m^). Une 



atténuation de 3 dB de part et d'autre de cette valeur maximale correspond à la "bande 
passante" du système. On a remarqué que cette bande passante s'élargissait quand 
l'amortissement augmentait. Si on appelle co x et 0) 2 les pulsations de part et d'autre de 

û) max qui correspondent à une atténuation de 3 dB de | Y (ûj)\ , on constate que 

20 log l0 = 20 log,„ = 20 log l0 -L = -10 log 10 2 = -3 dB 

1 max 1 max V ^ 

Cette relation permet d’écrire pour co i = OJ ] ou co i = Cû 2 

V 2 /y max = l/|y(^ )| => 42(2^ ) = -corf +( 2 Çco 0 co i ) 2 

qui en portant les deux membres au carré conduit à une équation du second degré en cor 

cof -2<w ( ,<»f(l-2^ 2 )+ry ( )(l-8^ 2 )=0, 



dont les racines sont 






De petites valeurs de Ç permettent d'écrire, en négligeant les termes en 



f 

1 + ^— 

2 v 



l + -^ + ^ 2 =1 + ^ 2 d’où tJ\ + Ç 2 =1 + — 
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ce qui conduit aux approximations suivantes 



03 

03 



2 

1 

2 

2 



= 03, 



03, 



(l-2f 3 )-2fVT7f î 



.(i-2f 2 )+2fVüi r r 



®»( l-2f)| 
®,;(i + 2f)i 



cù\ - 03; = AÇcol 



En considérant que O3 0 = 03 max et que 

03; - 03; = (o3 2 + 03, )(o3 2 — 03,) ~ 2 Û3 0 ( 03 2 - 03,) » 2(0 (û) 2 -03,), 

on obtient alors A 03 - 03 2 - 03, = 2 ty max Ç . La méthode consiste à faire varier la fréquence 
de part et d'autre du lieu de la valeur maximale pour noter les fréquences f, et f 2 qui 
correspondent à une atténuation de 3 dB, ce qui permet de calculer 



/ 2-/1 _ A/ _ A03 

(/ 2 +/ l )/ 2 /max ^max 



On définit aussi Q le facteur de surtension comme 



A 03 
03 



2 <r 



_i_ 

Q 




Figure 2.34 - Méthode de la bande passante à 3 dB 



Méthode d'identification du cercle ("circle-fitting") 

Au voisinage de 03 Q la fonction de transfert sur le plan de Nyquist est assimilable à un 
cercle (voir section 1.3). Ses paramètres (position et diamètre) sont identifiés à partir des 
valeurs mesurées par une procédure par moindres carrés et permettent d'obtenir ceux de la 
fonction de transfert H Fa ( 03 ) ou de l'admittance Y ( 03 ) : 03 Q , Ç et k (ou m ). 
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[ANNEXE] Informations pratiques sur les mesures vibratoires 



Calibrage des accéléromètres 

Système de calibration portable permettant de 
contrôler toute la chaîne de mesure (voir ci- 
contre). 

Des méthodes de calibration par comparaison 
avec un accéléromètre de référence (méthode dos 
à dos : ISO 5347-3) permettent des précisions de 
l’ordre de 1% à des fréquences fixes (80 Hz ou 
160 Hz) ou meilleures que 5% sur l’ensemble du 
spectre. Les méthodes de substitution utilisent un 
accéléromètre de référence intermédiaire ou un 
vibromètre laser. 




Frequency = 159.2 Hz 

co = 1 000 rad/sec 

Accélération = 10 ms -2 




Figure 2.35 — Système de calibration 
portable pour calibrer l’ensemble de 
la chaine de mesure. 



Standard Référencé 
Accelerometer 8305 
or Device Under Test 



Working Standard 
Accelerometer 4371 




Adaptor 

Plate 



Fixture 
WA 0567 




Figure 2.36 - Dispositifs pour V étalonnage des accéléromètres 



Montage des accéléromètres sur la structure 



■ par vis 

■ ciment exopy ou 
cyanoacrylate 

■ adhésif double face 

■ cire (pour des températures 
inférieures à 40°C 




Figure 2.37 - Modes de fixation des accéléromètres 
(Doc. Brüel & Kjaer) 
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Précautions d’emploi des accéléromètres 



Isolation électrique (rondelle de mica) de la 
structure pour éviter les boucles de masse 



Bruit triboéléctrique créé par les chocs du câble 
contre la structure 

Remède : fixer le câble à la structure par un 
ruban adhésif. 



Electrical 

(Prévention of ground loops) 




+ 40 
+ 30 

Les chocs violant comme les chutes sur les sols dB 
durs peuvent provoquer des accélérations de +20 

plusieurs centaines de g. L’endommagement +10 

est caractérisé par des altérations du spectre 
(d’après BK) 0 

- 10 

Figure 2.38 - Précautions à prendre pour 
l’ utilisation des accéléromètres (Doc. Briiel 
& Kjaer) 




Vibromètre laser 

La vibromètrie laser Doppler (LDV) est une technique de mesure sans contact des 
vibrations qui utilise l'effet Doppler. Le vibromètre permet la mesure sur des surfaces 
chaudes, molles, mobiles, miniatures, légères ou difficilement accessibles, sans les charger 
par la masse d'un accéléromètre. Les vibromètres utilisent généralement un laser Hélium- 
Néon (rouge visible A = 0.6328 pm) de faible puissance (classe 2) et peuvent être utilisés à 
des distances de quelques centimètres à plusieurs centaines de mètres. Ce sont des 
dispositifs interférométriques à deux faisceaux qui détectent la différence de phase entre 
une référence interne et le faisceau de mesure réfléchi sur la cible 3 . 

C’est donc la composante de la vitesse dans la direction du faisceau qui est mesurée. Les 
vibromètres délivrent souvent un signal analogique correspondant à cette vitesse vibratoire. 
Des systèmes comportant 3 têtes vont permettre de mesurer les trois composantes de la 
vitesse en un point de la surface de l'objet. 



pour des informations complémentaires sur ce matériel et son principe voir par exemple : 
www.polytec.de/polytec-com/l_vib/vib_uni_vib.html et www.bksv.com/pdf/Bpl864.pdf 
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Les vibrometers laser fonctionnent avec 
une grande variété d'états de surface mais il 
est important que la quantité de lumière qui 
retourne vers l'objectif soit suffisante. 
Quand la réflexion est spéculaire (surface 
polie), l'alignement doit être précis pour 
que le faisceau réfléchi revienne dans la 
direction de l'instrument. Au contraire, 
avec une surface diffuse, la lumière sera 
réfléchie dans toutes les directions mais 
parfois trop affaiblie pour les surfaces 
absorbantes. En cas de problèmes de 
mesure, il est possible d'utiliser une 
peinture ou un adhésif réfléchissant et 
diffusant (micro billes). 




Figure 2.39 - Vibromètre laser et analyseur de 
spectre {Doc. Polytec). 




Figure 2.40 - Un vibromètre à balayage est un système dont la tête comporte des miroirs pilotés 
par calculateur permettant d’effectuer automatiquement des mesures sur un maillage prédéfini 
(Doc Polytec). 
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LIBERTE 
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APPLICATION 
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4. INTRODUCTION A L'ANALYSE MODALE EXPERIMENTALE 

4.1 - Principe de l'analyse modale expérimentale 

4.2 - Le modèle modal 

4.3 - Acquisition des données expérimentales 
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ANNEXES 
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FONDAMENTAL 

1 - SYSTEME A DEUX DEGRES DE LIBERTE 



1.1 - Modèle non-dissipatif (non-amorti) 



h(t) 



x 2 {t) 



-> 






k 2 {x 2 - x,) 




AAAH 




- k x x x 














■ k 2 (x 2 - x x ) 



En sommant les forces qui s'exercent sur chaque masse, on obtient 

m x x x = -k x x x +k 2 {x 2 -Xj) 



m 2 x 2 = —k 2 (x 2 -x,) 



où encore sous la forme de deux équations différentielles du second ordre couplées 



m x x j + (k x + k 2 )x x - k 2 x 2 = 0 
m 2 x 2 + k 2 x 2 - k 2 x x = 0 



A ces équations correspondent 4 conditions initiales 

Xi (O) = Xjo , x 1 (0) = x 10 , x 2 (û) = x 20 , x 2 (û) = x 20 

Pour résoudre ce problème, on utilise une représentation matricielle qui présente l'avantage 
de pouvoir se généraliser à un nombre quelconque de degrés de liberté. 

Des vecteurs colonnes sont utilisés pour décrire déplacements, vitesses et accélérations 

x x {t) 



i{t) : 



x, (f) 

'2 



1 

X, {t) 



x{t)- 



X 2 (0 



x(t) : 



Xlit) 

'2 



1 

X, (t) 



La matrice des masses M 



m, 



0 



0 ni ~ 



permet d'écrire l'équation dynamique 



Mx 



m x 


0 " 






m x x x 


0 


m 2 _ 


1 

X: 

to 

1 




m 2 x 2 
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La matrice des raideurs permet d'exprimer les forces qui s'appliquent sur les masses 



K 



k ] + k 2 —k 2 
-k-, k , 



Kx = 


k l + k 2 —k 2 






(*1 +k 2 )x, -k 2 x 2 




~k 2 k 2 _ 






-k 2 x x +k 2 x 2 



Le système d'équations qui décrit le système à deux degrés de liberté peut donc s'écrire 
sous forme matricielle 

Mx + K x = 0 



Les matrices de masses et des raideurs sont symétriques 

M t =M et K t =K 

On suppose une solution harmonique complexe de la forme x(t) = X e jox avec X le 
vecteur des amplitudes indépendant de t (dans ce chapitre on omet volontairement de 
prendre explicitement la partie réelle : x(t) est donc une variable complexe). L'équation 
s'écrit 

(-CD 2 M + K)Xe ; “ =0 



Comme en général e jü * doit être différent de 0, l'équation précédente devient 

(-CD 2 M + k)x = 0 ou D X = 0 



D est nommée matrice dynamique. 



Une solution triviale correspond à D 1 DX = 0=>X = 0, (D 1 inverse de D), c'est à dire 
que = x 2 = 0 . La solution non-triviale correspond au cas où la matrice D est singulière 
(son déterminant est nul et elle n'a pas d'inverse). 



Rappel : si D = 



a b 


qlnrc — ci rf — lie pt ^ — 


" d 


-b 


c d 


dlUlo UCllI-/ / CLCC t/C CL , -, 

det(D) 


-c 


a 



La condition de singularité est appliquée à la matrice dynamique pour obtenir une solution 
indépendante de X : 



det(-&> 2 M + K)=0 



permet de déterminer l'inconnue co 2 



det(ü) = det 



CDm x +k l + k 2 

- k , 



- CD m 2 + k 2 



■ m 



x m 2 (D A - (»i^, +m 0 k 1 +m 2 k 2 )co 2 + k x k 2 =0 
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C’est une équation du second degré en cor qui fournit deux solutions co\ et (û \ , soit ± co l 
et ± (û 2 . En reportant cette solution dans l'équation matricielle, il est possible de 
déterminer deux vecteurs (vecteurs propres) X 

(-caf m+k)x 1 =0 
(-a> 2 M + k)x 2 =0 

Les vecteurs X, et X 2 ne sont connus qu'à un facteur multiplicateur près. La solution de 
x(t) est une combinaison linéaire des termes X x e ±ja '‘ et X 2 e ±ja>2 ' 

x(t) = [b i +B 2 e" M, )x i +(s 3 e j6>2 ' + B 4 e~ j<0l ‘)x 2 

où les coefficients B l à B 4 dépendent des conditions initiales x 0 = x(0) = [x lo x 20 ]' et 
x 0 = x(0) = [x 10 x 20 ] r . La solution peut aussi s'écrire sous la forme 

x(t) = A, sin(<yjf + ^i)Xi + A 2 sin(ft» 2 t + ^ 2 ) X 2 
où les amplitudes et les phases dépendent des conditions initiales. 

Chaque masse oscille aux deux fréquences û) l et co 2 appelées fréquences naturelles du 
système. Si les conditions initiales sont choisies telles que A 2 = 0 , alors les deux masses 
oscillent uniquement à (O x et leurs positions relatives dépendent du vecteur X , appelé 
déformée du premier mode. Réciproquement si les conditions initiales sont choisies pour 
que Aj = 0 . 

NB : les fréquences naturelles du système couplé à deux degrés de liberté ne 

correspondent pas à celle des systèmes isolés ou bloqués. 



1.2 - Exemple 

Soit le système précédent avec m, = 9, m 2 =1, k x = 24, k 2 = 3. L'équation 
caractéristique 

det(D) = ry 4 -6&r + 8 = (<w 2 -2 )(co 2 — ù)= 0 



conduit aux solutions Cüf = 2 et (û\ = 4 . 



Premier vecteur X 2 correspondant à car = (0{ = 2 



(-û)f m + k)x x 



"27-9(2) 


-3 ' 


*11 




" 0 " 


-3 


3 — (2) 


_* 12 J 




0 
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L'équation précédente fournit deux équations 9x n -3x n =0 et - 3x n + x l2 = 0 qui 
sont dépendantes et conduisent donc à la même solution 



x 



h 




donc le vecteur est défini à une constante multiplicative près 



X = 



1 /' 

73 

1 



Deuxième vecteur X , correspondant à or = CO 2 , = 4 

(-æ>2 m+k)x 2 



27-9(4) -3 " 

-3 3-(4) 



v 2\ 



22 . 



- 9x 2] — 3jc 22 = 0 soit x 2l = — — x 22 donc X 2 



Solution pour le mouvement des deux masses 



it) = 



x. 



w 



1 

X 2 {t) 



= A l sin (û) l t + </) l ) 



+ A 2 sin(&> 2 r + A) 



- 1 / 



avec û) l = V2 et û) 2 = 2 . Aj , A 2 , (j) x et (f) 2 sont déterminés à partir des conditions 



initiales. 



Solution correspondant à des vitesses initiales nulles : x 0 



1 

• 

o 

I 




1 

0 

1 


1 

*• 

to 

0 

1 







0 = (û x A 1 cos </) j x u + co 2 A 2 cos (j) 2 x 21 
0 = û) l A l cos A x 12 +co 2 A 2 cos (j) 2 x 22 



&>iA cos A 



= 0 



V*21 X 22 J 



71 71 

d'où cos A = 0 => A =-. De la même manière, on obtiendra (j) 2 = — . En prenant en 

compte les valeurs du déplacement initial, il est alors possible d’obtenir A, et A 2 



x 10 = Aj sin A x n + A 2 sin A, x 
x 20 = Aj sin (f ) j x p + A 2 sin (f) 2 x 



*io = A *n + A * 
X 



20 



21 

A *12 A *22 



soit 



^ _ * 10*22 * 20*21 
* 11*22 _ * 12*21 



et 



A, 



* 10*12 * 20*11 
* 21*12 _ * 22*11 



et finalement avec sin 



a+- 



K 



V 



= cos a , l'équation du mouvement dans ce cas particulier est 



‘(0 : 



M 



x 



*2 00 



— Aj cos (V2 r)-U cos(2 t ) 
Aj cos (72,)+' A 2 cos( 2 t) 
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Remarque : Si on veut mettre en évidence seulement le premier mode (à co l = V2 ), il faut 
que A 2 = 0 (c'est à dire x 10 x 12 =x 20 x n ), donc que les déplacements initiaux vérifient la 

relation — ^ . 

*20 *12 

Il en sera de même pour le deuxième mode (à û) 2 = 2 ), pour lequel A , = 0 conduit à 

*10 _ *21 
*20 *22 

Donc pour obtenir un mouvement correspondant à un mode unique, il faut que le 
déplacement initial soit proportionnel à la déformée modale du mode si les vitesses 
initiales sont nulles. 



1 .3 - Réponse forcée 

Des forces d'excitation F l (?) et F 2 (?) sont appliquées aux deux masses. L'équilibre des 
forces conduit au système d'équations 

M x(?) + K x(?) = F(?) 

Le second membre est le vecteur des forces F(?) = [F, (?) F 2 (/ )] 7 . Pour des excitations 
complexes F l (t)=F_ 1 e JÜ * et F 2 (?) = F 2 e J0 * (Fj et F 2 sont des nombres complexes dont 
les phases représentent celles de l'excitation). Au vecteur complexe F(?) = Fe J " Mf 
correspondra une réponse également complexe x(?) = \e i0M . Le système d'équations peut 
alors s'écrire 



D x = F 



avec la matrice dynamique (réelle dans le cas d'un système non-anorti) 



D 



■ ù) 2 m x + k x +k 2 —k 2 



- k 1 



— û)~m 1 + k 2 



La solution particulière s'écrit 



x = D 1 F 



où D 1 est l'inverse de la matrice dynamique D 

— co 2 m 2 + k 2 k 2 

k 2 — û) 2 m l +k x +k 2 

en dehors des fréquences naturelles (û [ et (û 2 où det(ü) = 0 les réponses complexes sont 
définies par 



D 



det(D) 
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2 - GENERALISATION AUX SYSTEMES A PLUSIEURS DEGRES 
DE LIBERTE 

Pour un système à n degrés de liberté, le vecteur déplacement devient un vecteur colonne 
«xl et les matrices de masse et de raideur sont des matrices carrées n x n . En considérant 
une masse j , l'équation d'équilibre s'écrit : 




(a) pour les masses j - 1 et j + 1 bloquées 



(b) pour les masses j - 1 et j + I en 
mouvement 



m j x j =-k j x j -k j+l x j 



■i=-k i U / - x , j ) - k (x - x /+1 ) 




avec x 0 = 0 , k n+l = 0 . L'ensemble de ces équations peut se mettre sous la forme d'un 
système linéaire 

M x + K x = 0 

avec 

M = diag(ra, , m 2 , m 3 , . . . , m n _ x , m n ) 



k 1 +k 2 -k 2 0 0 0 

— kr, k-, —k ^ '• ■ ■ 

0 -k 3 •. •. 0 : 

: o -k n _, o 

1 i ~k n _ i k n _,+k n -k 

0 0-0 -k m k n 



83 





III- Systèmes à plusieurs ddl 



Vibrations - Acoustique 1 



où M = diag(»z . ) matrice diagonale n x n pour un système à n masses. 

K matrice symétrique n x n des raideurs 

M t =M et K t =K (symétrie) 

Remarque : Dans le cas où la dernière masse est connectée par une raideur k n+l à une partie 
fixe, le dernier élément sur la diagonale de la matrice de raideur devient K nn = k n + k n+l . 

2.1 - Méthode de la base modale : réponse libre 

L'objectif est de ramener l'étude à celle d'un système à 1 DDL 

• en normalisant l'équation du mouvement par rapport à la masse, 

• en réalisant une transformation de coordonnées pour se placer dans la "base 
modale" où les équations du mouvement sont découplées. 

Soit le système à plusieurs degrés de liberté défini par 

Mx(f)+Kx(f) = 0 

dont on veut déterminer la réponse libre pour les conditions initiales: 
x(0) vecteur des déplacements initiaux à t = 0 , 
x(0) vecteur des vitesses initiales à î - 0 . 



Première étape : Normalisation par rapport à la masse 



V. v. 



La matrice des masses peut s'écrire M = M /2 M /2 où 



M /2 = 



m, 0 

0 ,/m 



0 0 



0 

0 

frn. 



= diag(^“). 



.-V. 



Il est possible d'exprimer M /2 = diag 
de variable) 



f \ 

1 



va/ w , y 



et de remplacer x par (premier changement 



.-V. 



x = M /2 q 



Le système d'équation du mouvement devient 



.-V. 



MM /2 q + KM /2 q = 0 



et en multipliant à gauche par M 



.-V. 



M MM q+M "^KM"^q = 0 
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On remarque que 

-Y -Y 

M /2 MM /2 = I : matrice unitaire 

_ 1/ _i/ ~ 

M 2 KM 2 = K : matrice symétrique normalisée de la raideur 
ce qui conduit à l'équation dynamique normalisée 

Iq+Kq = 0 

-Y i y 

Si x = M /2 q alors q = M 2 x et les nouvelles valeurs initiales sont: 

q(0) = M^x(0) et q(0) = M ! ^x(0) 



Deuxième étape : calcul des valeurs propres et vecteurs propres 

En cherchant une solution de la forme q(f) = \e l'équation du mouvement devient 

{-û) 2 I + K)\e j0X =0 

Après rejet de la solution triviale v = 0 et en posant A = (O 2 , on est ramené à un 
problème typique de recherche de valeurs propres : 

Kv = Av 

a) à chaque valeur propre A i correspond un vecteur propre v, , 

b) K est une matrice carrée n x n qui a n valeurs propres et n vecteurs propres, 

c) K est une matrice symétrique alors les A i sont réels et les v, sont composées 
d'éléments réels, 

d) les n valeurs propres sont obtenues en cherchant les solutions A de l'équation 
caractéristique 

det(AI-K)=0, 

e) les n vecteurs propres sont les solutions des n équations associées à chaque 
valeur propre A t 

U, i-kK=o 



f) les vecteurs propres forment un ensemble de vecteurs linéairement indépendants 
qui peuvent être normés. La norme du vecteur v est définie par 

i/ 



=V 



T 

y y = 



Z' 
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Pour normer le vecteur v , on recherche le scalaire a tel que u = a\ ait une 
norme unitaire, donc 

(crv T )(erv) = 1 => a 2 =—^— 

v v 



d'où 




g) Les vecteurs normés u, satisfont la relation d’orthogonalité 

c [0 i±j 



T 

u u 



1 



1 = J 



et peuvent définir une matrice orthogonale (les vecteurs u ; sont dits 
orthonormaux) 

P = K u 2 ••• uj 

telle que P T P = I et P T K P = diag(A f ) . 

Conséquence: P _1 =P T (voir Annexe) 



Troisième étape : projection dans la base modale 

La matrice orthogonale P est utilisée pour effectuer un nouveau changement de variable : 
on pose 

q(f)=Pr(f) 

et on multiplie à gauche par P T 



P T Pr(() + P T K Pr(f) = 0 

En tenant compte des propriétés énoncées plus haut et en considérant A t = co 2 , cette 
transformation conduit à l'expression suivante 

I r(?) + diag [co 2 ) r(t)=0 



qui représente un système de n équations indépendantes ( découplées ) : 

q ( t ) + ryj 2 r x ( t ) = 0 

r i {t) + (0 2 r i {t)= 0 

f„ (t) + ûj 2 r n (f) = 0 



86 






Vibrations - Acoustique 1 



I I I - Systèmes à plusieurs ddl 



En exprimant les conditions initiales dans la base modale 

q(0 = Pr(f) => r(t) = P _1 q(t) = P T q(t) 
d'où r(o) = P T q(û) et r(û) = P T q(0) 

la solution pour r i ( î ) est celle du système à un degré de liberté 



M V^oHo . f , , + 0 )^ o 

r ; \t) = sin co ( t + arctan — — 



( r 0 et r i0 sont les éléments des vecteurs r(û) et r(o) ) 

Quatrième étape : calcul des déplacements des masses ( transformation inverse ) 

En utilisant les équations précédentes, il est possible de calculer le vecteur des 

déplacements x à partir des solutions dans la base modale r . Puisque x(f) = M q(7) et 

q{t)=Pr{t) 

x = M ^q = M ^ 2 Pr = Sr 

La matrice S = M 72 P où P = [uj u 2 • • • u (I ] est donc constituée de vecteurs colonnes 

liU'X/Jrrh] 



v*i )\ / 

tels que q».=M"^u ( . = _ 

<ï>, est le vecteur de la déformée modale du mode i, c'est à dire le vecteur des 
déplacements unitaires (<!*)• de chaque masse j . 

Les déplacements peuvent s'écrire 

K 

x(;)=S r(/)- V/;(i)'t> , avec <t> = {i 

ï=l ; 

soit pour une masse j 



n ri 

i W =X r i (0 (* / )., =Z A- sin^f + (pi ) 
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avec A, 






2.-2 
iO + r i0 



CD, 



et (Ç; = arctan 



<°i r i0 



Résumé de la démarche 

Première étape : Normalisation par rapport à la masse et premier changement de variable : 

x — > q 

Deuxième étape : recherche des valeurs propres A j = cof et des vecteurs propres v ( . . 

Les vecteurs normalisés u ; obtenus à partir des vecteurs v i permettent de construire une 
matrice orthogonale P 



Troisième étape : projection dans la base modale 

La matrice P est utilisée pour effectuer un changement de base qui conduit au deuxième 
changement de variable : 

(1 > r 

Le système d'équations obtenu correspond à n équations différentielles indépendantes de 
systèmes à un degré de liberté dont les solutions r ( t ) peuvent être facilement calculées en 

employant les conditions initiales ( x 0 — > q 0 — >r 0 ,x 0 — >q 0 — > r 0 ). 



Quatrième étape : calcul des déplacements dans la base physique ( transformation inverse) 
Pour cela le changement de variable inverse est effectué : 

r — > x 

Le déplacement x ■ (t) de chaque masse est une combinaison linéaire des solutions r (t) de 
la base modale : 



,to=Er ( (f)(40,=£ 



i = 1 



i=l 






CD: 



+ 



f 2 

'iO 



r 



-sin 



CD, + arctan - 



ù), r 



v 



1 0 



«0 ) 




2.2 - Systèmes avec amortissement visqueux 

Les systèmes réels sont amortis mais on ne connaît pas bien dans la plupart des cas le 
modèle d'amortissement. Souvent, le modèle d'amortissement visqueux est utilisé pour des 
raisons de simplicité. La méthode consiste à considérer un taux d'amortissement modal 
< 1 inclus dans les équations découplées de la base modale 

r{t)+2C i Cü i r i {t) + C0fr{t)= 0 
La solution est (voir système à 1 ddl) 
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r t ( t ) = A,. exp(- Ç^t) sin(û) a ,t + à, ) 



avec 



\ 



(r i0 +Ci a^io ) 2 + (no® ai) 



,,2 
<2 / 



et (f) i = arctan 






A + 



Le taux d'amortissement modal devra être estimé expérimentalement, soit être identifié 

à partir de la matrice des coefficients d'amortissement qui entre dans le système d'équations 
du mouvement : 



Mx(t)+ Cx(f)+ Kx(f) = 0 



*i (0 



x 2 (t) 



-> 






A 



^vw4 

I — 1- 



m, 



KV\aH 

S- 



m , 



Par exemple, pour un système à 2 ddl 



m l 


0 ' 


X + 


Cl +C 2 


-c 2 l 


X + 


k x + k 2 -k 2 


0 


m 2_ 




_ -c 2 


r 

c 2 J 




_ -k 2 k 2 _ 



Cette équation pose problème car la décomposition modale ne peut pas être réalisée par la 
procédure décrite précédemment car les coefficients d'amortissement créent un couplage 
supplémentaire entre les équations du mouvement. Plusieurs méthodes numériques sont 
utilisables pour résoudre ce problème. Cependant, dans bon nombre de cas, il sera possible 
d'approcher la matrice C par une combinaison linéaire de matrices de masse et de raideur 

C = «M + /?K, a et j3 constants. 

Cette forme d'amortissement est appelé amortissement proportionnel 

M x(t) + [crM + /?K]x(/)+K x(?) = 0 

En posant x(t) = M /2 q(t) et en multipliant à gauche par M , il vient 

Iq(t)+[aI + ^Kjq(t)+Kq(t) = 0 
puis en posant q(t) = Pr(/) et en multipliant à gauche par P 1 

I r(t) + [a I + p diag(&> 2 )]r(f) + diag(&> 2 ) r (t) = 0 
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Cette dernière équation ne comporte que des matrices diagonales donc correspondant à n 
équations modales découplées 

ï, ( t)+ [oc + fîojf ] r. (?) + cof /• (?) = 0 

En rapprochant cette équation de celle utilisant le taux d'amortissement modal , on 
obtient l’équivalence 2 Ç i (û i = a + j3cof, d’où 

„ a fico . 

Ci = — + ! — L . 

2 co. 2 



2.3 - Réponses forcées 

L'équation se met sous la forme 

Mx(?)+ Cx(?)+ Kx(?) = F(?) 

F(?) est le vecteur des forces appliquées à chaque masse 

F(?) = [Fj(?) F 2 (t) - F n (?)] 7 

En suivant toujours la même démarche : on pose x(?) = M 22 q (/ ) et on multiplie à gauche 

_i/ 

par M /2 

I q(?) + C q(?) + K q(?) = M “ F(?) 

avec C = M 72 CM 72 , puis on pose q(?) = Pr(?) et on multiplie à gauche par P 1 

If(?)+diag(2^?y,.)f(?) + diag(?y i 2 )r(?) = P T M” 2 " 2 F(?) 
où est obtenu par exemple par la méthode de l’amortissement proportionnel. 

f(?) = P A M 22 F(?) est le vecteur des forces modales (généralisées), dont les éléments 
/, (?) sont des combinaisons linéaires des forces F i (?) . Finalement, l'équation modale 
découplée est de la forme 

Ï, (0 + 2 Ci co, r, (r) + CO , 2 /• (?) = f; (?) 

dont la solution complète est celle du système à 1 ddl de déplacement r (?) 

r, (/) = A exp(- C;CO i t)sin(co a ,! + </),)+ f, (?) * h, (?) 

1 1 

= A exp(- C;CO;t) sin(© fl ,? + </>,) + exp(- Ç^t ) f /, sin co rn (t-r)dr 

c o„ { 
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Les coefficients A f et (f) i de la solution générale sont déterminés en utilisant les conditions 
initiales. 

En régime stationnaire (harmonique permanent), la solution particulière pour 
F j ( t ) = F j0 cos (tôt + (p j ) est dans la base modale 

r (t) = = cos O) t - arctan + 9 i 

yj(a)f-(0 2 ) 2 +{2 GXO.Q 2 \ ~ ùy ) 

avec f 0 = P T M"^ F 0 = S T F 0 où S T = (m~ 1/2 p) r = P T M" 1/2 . 

f 0 est le vecteur des forces modales généralisées complexes dont chaque élément est 
/ = f i0 e Jd ‘ ■ Il est calculé à partir du vecteur F 0 des forces complexes appliquées à 

chaque nœud (masse) du système et dont chaque élément s'écrit F_ j0 = F j0 e iç>i . 
L'utilisation d'un vecteur complexe est nécessaire pour représenter les déphasages (p j 
différents pour chaque force. 

Le vecteur déplacement s'écrit comme précédemment par 

x(t)=S r(f) = £#;■(*)<&,■ 

Ï=1 

NB : En utilisant la fonction de transfert modale complexe L {co) = -, — rp , 

[cor -m-)+2j(Oû) i Ç i 

la solution dans la base modale s'exprimer simplement par r i (t) = Re{/ L ( . (oj) e ><0 ' }. 
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APPLICATION 



3 - ABSORBEUR DYNAMIQUE (système à 2 DDL) 

C'est un système employé pour réduire les vibrations d'un premier système assimilable en 
basse fréquence à un système à un ddl (une machine montée sur une suspension élastique, 
par exemple). Le problème vibratoire se pose en général à la pulsation de résonance du 

premier système Çl l =^Jrn~ l (au cours d’une montée en vitesse). Pour réduire 
l'amplitude des vibrations, une masse m 2 est suspendue élastiquement par une raideur k 2 
au premier système (la machine) constituant ainsi un système à deux ddl. 




Problème quand û) = Q. 1 Ajout de l'absorbeur dynamique 



A partir des relations précédemment établies pour le système à 2 ddl, avec F 2 = 0 , les 
amplitudes des deux masses sont 




En posant Q.f = 



(pulsations naturelles des systèmes lddl isolés) 




Cette équation montre que quand cû correspond à la fréquence naturelle de l'absorbeur 

nr 

seul 0. 2 = — , l'amplitude des vibrations de la masse 1 s'annule : x, = 0. On accorde 

V m 2 

donc ü. 2 sur la fréquence dont on veut réduire l'amplitude. 



Souvent, la pulsation initiale à réduire est co = £l l = —, l'absorbeur est accordé sur 

"\ m, 



cette valeur 
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Utilisation d'un absorbeur dynamique pour réduire l'amplitude des vibrations 
d'un système à sa pulsation propre 



4 - ANALYSE MODALE EXPERIMENTALE 



4.1 - Principe de l'analyse modale expérimentale 

La représentation modale permet de synthétiser le comportement dynamique des structures 
comme la somme de plusieurs systèmes à 1 ddl. L'intérêt de cette représentation trouve sa 
justification dans le fait que les amplitudes de vibration deviennent très importantes autour 
des fréquences de résonance de ces systèmes. 



Ainsi, le test vibratoire d'un système mécanique aura souvent pour objectif 

■ de connaître ses fréquences propres 

■ de déterminer ses autres paramètres (amortissements, vecteurs propres) pour 
obtenir un modèle analytique théorique complet qui pourra ensuite être 
employé en simulation pour connaître son comportement dynamique avec des 
excitations particulières. 

C'est l'objet de l'analyse modale expérimentale. Elle est basée sur les approximations 
suivantes : 

■ le système est linéaire (dans la gamme des amplitudes étudiées), 

■ le système continu peut se représenter par un système discret où les paramètres 
sont exprimés pour chaque nœud du maillage (nombre de ddl = nombre de 
mailles), 

■ l'amortissement est supposé proportionnel. 
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4.2 - Le modèle modal 

4.2. 1 - Réponse d’un système à N ddl 

Soit un système à N ddl dont le modèle peut se décrire par la relation matricielle suivante 
(Rappel): 

Mx(/)+Cx(f)+Kx(t) = F(t) 

avec le vecteur des déplacements x(/) = [x l , ■ ■ ■ , x N ]' aux différents points du maillage et 
des forces F(/) = [F l , ■ ■ ■ , F N ] T . En notation complexe 

[- &rM + j(0 C + k] x = F 

(x et F sont complexes avec une dépendance temporelle en e j0> ‘ ) 




La recherche des valeurs propres de la matrice de raideur normalisée par la matrice de 
masse M = diag(/;z y ) conduit aux fréquences propres et aux vecteurs propres : 



K = M 



km"^ 



valeurs propres 
vecteurs propres normés 

{Uj.-.U,.,-,!!*} 
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Les vecteurs propres normés constituent une matrice orthogonale P = [u x • • • u, , • • • n N ] qui 
permet de construire la matrice de passage entre la base physique et la base modale : 



r = S 1 x 







x = Sr 



Le système d’équations découplées dans la base modale 

I r(t)+ diag(2£.®,. )r(t)+diag(<sr )r(/) = f (t) 
qui peut aussi s’écrire en notation complexe 

[- co 2 l + diag (j2Ç i o) i û))+ diag (cor )]r = f 
avec le vecteur des forces modales généralisées 

f = P T M % = S T F 



Les solutions de ces N équations sont 



avec 



M = /, M A (oj) 



r (oj) déplacement modal complexe pour le mode i, 

f i (a>) force modale généralisée complexe appliquée au mode i, 

L i (oj) fonction de transfert complexe pour le mode i 



L,.(ty) = 



cor 



1 

(O 1 + Jjl^ t (O t (O 



Les déplacements aux points du maillage s’obtiennent par 



x(&>) = S r(<») 



(<y) = X r t (oj) <î> ,■ = X fi (ûj)L i (oj) <î> 

1=1 1=1 



<ï* ; est le vecteur qui correspond à la déformée modale de chaque point du maillage pour le 
mode i . Chacun des N vecteurs sont les colonnes de la matrice S = [<D,, 
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4.2.2 - Expression des fonctions de transfert 

En remplaçant le vecteur des forces modales généralisées par son expression f = S T F 



x = S diag(L ; ) S T F = ^ — 

tï CQ- 



co 2 + j2Ç i co i co 



F 



Si on veut exprimer la réponse x m (ûj) à une excitation F n (co) en un seul point n du 
maillage, le vecteur force F = [0, • • • ,0, F n (co),0, ■ ■ ■ ,0] est utilisé dans l’expression 
précédente 



»=I 



(4»,.) (0»T ) 

V i / m V i / n 

co 2 -œ 2 + jlÇ'OOtû) 



f.M 



En définissant l’admittance Y mn {(û) comme la fonction de transfert x m {(û)/ F n (oj), on 
constate que la matrice des admittances est 



Y Y 

-Ml 1 12 

Y Y 

1 21 1 22 



IN 



<ï> 4» 1 

i i 



y - 

tt O) 2 - co 2 + j2( i û) i û) 



dont les éléments sont des fonctions de transfert entre deux points 



r„» 



■*„(<») 

F» 



g (»,)■ (*1). 

m ty, 2 - « 2 + j2Ç i m i co 



La symétrie de la matrice Y est une propriété connue comme le principe de réciprocité 



Y = Y t 



Y (co) = Y (co) 

mn V / nm V / 



x M) 

F.M F.M 



F m M 





96 



Vibrations - Acoustique 1 



III- Systèmes à plusieurs ddl 



4.3 - Acquisition des données expérimentales 

L’analyse modale expérimentale est une méthode inverse : à partir des réponses en chaque 
point du maillage à des excitations (forces) connues, nous cherchons à déterminer les 
caractéristiques du modèle à N ddl. 

2.3.1 - Conditions d’essai 

Le mode de fixation de la structure est très important. Le plus souvent la structure est : 

■ suspendue (libre), 

■ encastrée, 

■ en condition de fonctionnement. 

Selon la configuration choisie, les modes risquent d’être différents (en particulier les 
premiers). Les fixations libres ou encastrées permettront des comparaisons avec des calculs 
par éléments finis. 




Figure 3.1 — Excitation d’une structure d’avion par des pots vibrants pour déterminer 
des fonctions de transfert expérimen tales (Doc. ONERA) 
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4.3.2 - Mesure des fonctions de transfert entre les points du maillage 

On trouvera différentes formes en fonction de la réponse (déplacement, vitesse ou 
accélération) 



Accélérance 



An M : 



M 



F n M 



Mobilité 






F, (®) 



CO 



Admittance 



y (~) AJ®) 






ty 



ÛT 



avec l’accélération a m (co) et la vitesse vibratoire v m (co) = a m ( CO )/ jco . 



La fonction de transfert mesurée (FRF : Fonction de Réponse en Fréquence) dépend du 
type de capteur utilisé pour mesurer la réponse de la structure (accéléromètre — > a m (co), 

vibromètre laser —> v m (co)). Pour généraliser, nous désignons par H nm {co) la fonction de 
transfert mesurée entre un capteur placé au point m et une force excitant le point n : 
H mn {co)e {A mn (co\ M mn {co), F m „ (ty)}.Très souvent les logiciels d’analyse modale 
travaillent sur des admittances Y mn {co) . En fonction du capteur employé pour mesurer la 
fonction de transfert, elle devra être divisée par jco (vibromètre laser) ou par - co 1 
(accéléromètre). 





D Accéléromètre 
n Capteur de force 




Figure 3.2 — Les deux modes d’excitation couramment utilisés pour obtenir les FRF 
entre les points du maillage d’une structure : le pot vibrant (gauche), le marteau 
d’impact (droite). 



Principalement, deux types d’excitations sont employées (voir Chapitre II) : 

a) le pot vibrant 

Il permet d’exciter d’injecter d’avantage d’énergie, donc d’exciter des structures 
plus lourdes. En contre partie, son installation est délicate et il ajoute localement 
de la masse à la structure. Le signal d’excitation peut être un bruit blanc, un 
bruit pseudo-aléatoire synchronisé sur la fenêtre d’acquisition, un balayage de 
signal harmonique, un chirp, etc. 
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b) le marteau d’impact 

Il est facile à mettre en œuvre mais réclame une certaine expérience. Bien 
adapté aux structures légères, il ne les charge pas et le spectre d’excitation peut 
être ajusté par le choix d’un embout. Pour les structures très résonantes, il faut 
s’assurer que le signal de réponse présente une ‘extinction’ suffisante pour ne 
pas être tronqué à l’extrémité de la fenêtre d’acquisition, en lui donnant une 
dimension adéquate ou/et en appliquant une fonction exponentielle. Cette 
dernière procédure étant équivalente à l’ajout d’un amortissement 
supplémentaire, une correction devra être apportée par le logiciel, après 
détermination des paramètres modaux. 



4.4 - Extraction des paramètres modaux 

Il existe plusieurs méthodes pour F extractions des paramètres modaux à partir des mesures 
sur la structure dont les deux grandes catégories sont : 

■ les méthodes fréquentielles 

Des fonctions de transfert idéales de systèmes à un ou plusieurs ddl sont ajustées 
sur les courbes expérimentales par des méthodes de moindres carrés ( méthode de 
lissage des fonctions de transfert) 

■ les méthodes temporelles 

Ces méthodes sont basées sur la représentation du comportement à partir des 
réponses impulsionnelles modales. Les réponses sont obtenues en multi-excitation 
et les réponses impulsionnelles souvent calculées par transformée de Fourier 
inverse des fonctions de transfert ( méthodes d’ibrahim et des exponentielles 
complexes). 




Figure 3.3 - Fonctions de transfert obtenues expérimentalement entre les différents 
points du maillage 



D.J. Ewins, Modal Testing : Theory , Practice, and Application (2 nd édition), Research Studies Press, 2000. 



2 J. Piranda, Analyse modale expérimentale, volume « Mesures et contrôle », Techniques de l'Ingénieur, 
2001 . 
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